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Abstract 



In pursuit of the origin of the masses of particles in the Standard Model, the Brout- 
Englert-Higgs Mechanism predicts the existence of a scalar boson, whose discovery has 
been the main goal of the Large Hadron Collider. However, introducing this scalar boson 
does not provide the answer to all inadequacies of the Standard Model. Many different 
Beyond Standard Model theories have been proposed in order to explain these anomalies. 
Multi-Higgs-doublet models are among these models. 

We restrict ourselves to the scalar sector of multi-Higgs-doublet-models. The huge 
number of free parameters in the scalar potential in these models, makes it impossible to 
study the most general case for any N. In order to reduce the number of free parameters, 
one could impose symmetries on the scalar potential. Therefore, it is important to 
explore which symmetries can be implemented in the scalar sector, how these symmetries 
are broken, how they could be encoded in the Yukawa sector and what the resulting 
properties of the fermions are. 

Classifying these symmetries in the scalar sector is the main focus of this thesis. We 
have found certain symmetries that are always broken in models with more than two 
doublets, which we name "frustrated symmetries". Examples of such symmetries are 
presented in 3HDM, and one particular symmetry, the octahedral symmetry, is studied 
further. This symmetry seems to be the largest realizable discrete symmetry that can 
be imposed on the scalar sector in 3HDM, and interestingly results in a 2HDM-like mass 
spectrum. 

In the attempt towards the classification of possible symmetries in the scalar sector 
of the NHDM, we find that these symmetry groups are either subgroups of the maximal 
torus, or certain finite Abelian groups which are not subgroups of maximal tori. For the 
subgroups of the maximal torus, we present an algorithmic strategy that gives the full 
list of possible realizable Abelian symmetries for any given N. We extend this strategy 
to include Abelian antiunitary symmetries (with generalized CP transformations) in 
NHDM. 

We also show that multi-Higgs-doublet models can naturally accommodate scalar 
dark matter candidates protected by the group Z p , since these groups are realizable in 
NHDM. These models do not require any significant fine-tuning and can lead to a variety 
of forms of microscopic dynamics among the dark matter candidates. 

The results of this thesis are published in [I3J HHl [63] . 
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Chapter 1 
Introduction 



This is the most exciting time to be a physicist; With the recent and upcoming results 
of the Large Hadron Collider (LHC) physicist are gaining better knowledge on the prop- 
erties of the basic constituents of matter, and the fundamental interactions in nature. 
The current best-suited-with-experiment framework describing high energy phenomena 
is the Standard Model (SM) of electroweak and strong interactions. 

The first blocks towards the construction of the Standard Model were put down by 
Glashow in 1960, where he combined the [/(l)-symmetric theory of electromagnetic in- 
teractions and the 5't/(2)-symmetric theory of weak interactions into an SU(2) x U(l) 
symmetric electroweak theory pGj. Later on, in 1967 Weinberg [2] and Salam [3] incor- 
porated the Brout-Englert-Higgs mechanism [4 J [5J [6 J into Glashow's electroweak theory, 
and gave the model its modern form. And finally, in 1970's the S'?7(3)-symmetric strong 
interactions were described by quantum chromodynamics (QCD), shaping up the Stan- 
dard Model. The symmetry group of the Standard Model is denoted by SU(3)c x 

SU(2) L X U(l)y. 

Despite its great success in explaining all available data, the Standard Model has 
several inadequacies. It requires the existence of a massive scalar particle, the Brout- 
Englert-Higgs boson, whose discovery has been the main challenge of the LHC. Moreover, 
it does not give an explanation for the fermion mass spectrum, in particular the small 
neutrino masses (and their oscillations). The theory does not contain any viable dark 
matter particle that possesses all of the required properties deduced from observational 
cosmology. Also, this model does not include gravitatiorQ This theory also falls short 
in explaining the strong CP problem, and matter-antimatter asymmetry. 

These deficiencies led physicists to explore Beyond the Standard Model theories 
(BSM). A rich spectrum of theories has been proposed, such as supersymmetric ex- 

1 This is partly because at the level of elementary particles the gravitational force is many orders of 
magnitude weaker and can be ignored. The quantum effects of gravity become important in describing 
particle interactions at the Planck scale. This refers to either a very large energy scale (1.22 x 10 19 
GeV) or a very tiny size scale (1.616 x 10 35 meters). 
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tensions, string theory, extra dimension theories and grand unified theories, etc. The 
question of which theory is the one chosen by the Nature, can only be answered through 
experiments. 

In this Chapter we study the Standard Model Lagrangian and present the particle 
content of the SM in Section [LT Section L2 presents the Brout-Englert-Higgs mecha- 
nism through which the gauge bosons acquire mass. We end this Chapter with motiva- 
tions for introducing more than one doublet, and give the example of a two-Higgs-doublet 
model. 



1.1 The Standard Model of particle physics 

The particle content of the Standard Model is shown in Figure ( |1.1[ ). The three families 
of fermions (with spin 1/2), which obey the Fermi-Dirac statistics, form matter. Each 
family consists of two quarks and two Leptons. These families are a perfect replica of 
each other, where the heavier families eventually decay into the first family particles 
which are stable. Each fermion is associated with an anti-particle with the same mass 
and decay width but opposite charges. 

The bosons (with spin 1), which obey the Bose-Einstein statistics, mediate the funda- 
mental interactions. The photon 7 is exchanged in electromagnetic interactions between 
charged particles. The gauge bosons W + , W~ and Z° mediate the weak interactions of 
radioactive decay (of all fermions). The eight gluons are responsible for strong interac- 



tions inside the atomic nuclei (among quarks). The Feynman diagrams in Figure (1.2) 
illustrate these three fundamental interactions. 

There is also the two theoretical bosons, the graviton (with spin 2) which is re- 
sponsible for the gravitational interactions, and the Higgs boson (with spin 0), which is 
responsible for the masses of the SM particles. 



1.1.1 The Standard Model Lagrangian 

The Standard Model is based on a very powerful principle, the gauge symmetry; the 
Lagrangian remains unchanged under gauge transformations of an internal symmetry 
group, SU(3) C x SU(2) L x U(1) Y . 

The Standard Model consists of the electroweak theory describing the electroweak 
interactions which is based on a SU(2)l x U(1)y symmetry, and the theory of Quantum 
Chromodynamics (QCD) describing strong interactions with an SU(3)c symmetry. 

In this Chapter we only study the QED and electroweak Lagrangian. 

QED Lagrangian 

The Lagrangian of QED theory involves three parts; 

LqeD — Lem + Lf)i rac + Li nt (1-1) 
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fermions 



quark up 
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tau 
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7 
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z° 

o 
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±1 

80.4 GeV 



Higgs 

H 



> 114 GeV 
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Figure 1.1: The elementary particles of the Standard Model, and their characteristics 
0- 



The electromagnetic field Lagrangian describes the propagation of free photons; 

Lem = —-^F i ivF ilv 

with = d v A a - d a A u 



(1.2) 



where is the four-vector electromagnetic potential. The Dirac Lagrangian has the 
form: 

L D irac = iip^d^ip - m^ip (1.3) 

And the Lagrangian for interaction of a particle with charge q and the electromagnetic 
field is described with: 

L^t = -q^i>A a (1.4) 
This Lagrangian respects the U(l) gauge symmetry; 



if) -> e~ iqa{x) i) 



Au -> Au + -dua(x) 
q 



;i.5) 
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a 




Figure 1.2: The three fundamental interactions of the Standard Model [7]; (a) The 
electromagnetic interaction with the exchange of a photon, (b) the weak interaction 
with the exchange of a W boson, (c) the strong interaction with the exchange of gluons 

which is equivalent to replacing the normal derivative by the covariant derivative D^; 

d a — > , D ll = d ll + iqA ll (1.6) 

This replacement preserves the local symmetry. The electric charge q acts as a generator 
for the local group transformation U(x) = e~ iga ^ x \ According to Noether's theorem, 
the continuous symmetry results in a conserved quantity that is the electric charge. The 
massless gauge boson associated with the local U(l) symmetry is the photon with couples 
to particles with the electric charge. So, QED is an Abelian theory with the symmetry 
group U(l) and one gauge field, the electromagnetic field. 



The electroweak Lagrangian 

The electroweak Lagrangian is an SU(2)l x U(l)y gauge invariant theory. There are 
three gauge bosons associated with the SU(2); W^, W^. There is an additional 
gauge field associated with the U(l) symmetry B^. The electroweak Lagrangian has the 
following form; 

L Ll e pt Qn ~\- Lg aU g £ 

Llepton = i'lpRl^iP^R + i'lpjffip^L 

L 9 au ge = -\fa„r -\Tr{F, u Fn (1.7) 
where the strength tensors of W and B fields are; 

ffiu = d^B u — d u B u 
daW l v - d v W\ - g w e lmn W™W? 
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and the covariant derivative preserves the local gauge invariance; 

D„ = d„ + %1 ^B, + i 9w ^- ■ (1.9) 
Here 7^ are the gamma matrices (Dirac matrices), and Oi are the Pauli matrices 

(Ti = ( ? J ) , 02 = ( °, * J , ^3 = ( n ° J (1.10) 



x 1 y z y z j ' d y -1 y 

and are the SU(2) and C/(l) gauge couplings respectively and e lmn is the totally 
antisymmetric symbol. 

Note that the electron and its neutrino field are combined together into a doublet for 
the left-handed components, and a singlet for the right-handed component; 



^ L= [e L J' ^R = e R (1-11) 

since the mass of the neutrino is in SM and there is only a left-handed component of 
the neutrino; 

1-75' 



: = Ve (1.12) 

The left-handed and right-handed electron fields, €l and e R , are related to the electron 
in the usual way; 

e L = e, e R = [ e (1.13) 



The three types of charges in the electroweak theory are the electric charge Q, the 
weak isospin /, and the weak hypercharge Y. These charges are related by the Gell- 
Mann— Nishijima relation; 

Q = I 3 + ^ (1-14) 

where I 3 is the third component of the weak isospin. The neutrino is assigned an isospin 
of I 3 = +1/2. A left-handed electron has If = —1/2, and a right-handed electron field 
has I R — 0. For a left-handed spinor Y — — 1, while for a right-handed spinor Y = — 2. 
So, the overall charges are: 



Neutrino 


T3 - ,1 
1 L — + 2 




-1 


Ql 


= 


Left-handed electron 


= 3 = -| 


n = 


-1 


Ql 


= -1 


Right-handed electron 


: J3 = Q 


n = 


-2 


Ql 


= -1 



(1.15) 

The four gauge fields W^, W^, W 3 and correspond to the weak isospin charge / and 
weak hypercharge Y, respectively. Clearly a particle with a given type of charge, interacts 
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with the field associated with that charge, and the value of the charge determines the 
strength of that interaction. 

The left-handed fermionic fields with weak isospin 7=1/2 (I 3 = ±1/2) and weak 
hypercharge Y{Q L ) = 1/3 can be written as: 




(1.16) 



For the fermions of right helicity, the weak isospin is zero (I3 = 0) which corresponds to 
singlets of fermions: 

Ur 2 ' 3 = u R ,c R ,t R , D R ' 2 ' 3 = d R ,s R ,b R (1.17) 



1.2 Spontaneous symmetry breaking 

The local gauge symmetry causes the W and Z boson masses in the electroweak La- 
grangian to be zero, which is far from the experimentally observed values, = 80GeV 
and Mz = 91GeV. Introducing a mass term, by hand, to the Lagrangian is not the 
solution, since it breaks the gauge symmetry. The mechanism that insures massive W 
and Z bosons and massless photons while preserving the gauge invariance, is called the 
Brout-Englert-Higgs mechanism, and is presented in this Section. 

This mechanism also explains the masses of the quarks and leptons through Yukawa 
interactions with the scalar fields and their conjugates. 



1.2.1 Abelian Higgs model 

Let's start with an Abelian Higgs model, where we introduce spontaneous symmetry 
breaking in a model with a (7(1) symmetry, say electromagnetism theory. 
Recall that the [/(l)-symmetric electromagnetism Lagrangian; 

Lem = —-^FjxvF^ 

with F llv = d v A ll -d li A v (1.18) 

is invariant under the transformation A^ — > ^d^a(x). Let's add a mass term to the 
Lagrangian: 

L = ~F^ + l -m 2 A^ (1.19) 

The mass term clearly violates the gauge symmetry. However, we can extend the model 
by adding a single complex scalar field; 

0= -L(0 1+ ^ 2 ) (1.20) 
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with real fields <j>\ and <p2- The Lagrangian describing this new field's kinetic and potential 
term has the form 



F F^ u + ID 



with V(0) = [i 2 
and D^ = + iqA^ 



2 " ^(0) 

2 + A(|0| 2 ) 2 



:i.2i^ 



where V(4>) is the most general [/(l)-symmetric renormalizable potential. For A > 0, 
this theory has two different forms: 



/i 2 > 

In this case the potential has a shape as shown in Figure (1.3), and it preserves 
the symmetry of the Lagrangian. The vacuum appears at = 0. This Lagrangian 
describes a massless photon and a charged scalar 0, with mass \x. 









4>' 



Figure 1.3: The scalar potential with /x 2 > 



fi 2 <0 

The more interesting potential in this case has the shape shown in Figure ( 1.4 ) . The 

vacuum appears at (0) = \J — which breaks the symmetry. It is conventional 
to choose the vacuum to lie along the real axis of 0. 




Figure 1.4: The scalar potential with /z 2 < 
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Goldstone bosons 

It is common to parametrize the field in terms of two real fields x an d H, 

<t>=-L e ^{v + H) (1.22) 
The Lagrangian in terms of x an d H, will have the form; 

+^H 2 + l -q 2 v 2 A,A» 

—qvA^d^x + higher order terms (1-23) 

This Lagrangian describes a photon with mass M4 = qv, a scalar filed with mass 
a/— 2/i 2 > 0, and a scalar field with mass zero. The confusing mixed term —qvA^d^x 
can be removed by making a gauge transformation: 

Afj, > A^ — -d, X (1-24) 
qv 

This transformation removes the field x from the Lagrangian. The field x is called a 
Nambu-Goldstone boson P [9], which has been "eaten" to give mass to the photon. The 
field H is called a Higgs boson field. 

Let's count the number of degrees of freedom (dof). Before the spontaneous symmetry 
breaking we had a massless photon (2 dof) and a complex scalar field (2 dof), making 
the total of 4 degrees of freedom. After the spontaneous symmetry breaking we have a 
massive photon (3 dof) and a real scalar H (1 dof), making the same total number of 
degrees of freedom. 



1.2.2 Non-Abelian Higgs model (Brout-Englert-Higgs mecha- 
nism) 

Here we present the spontaneous symmetry breaking method which forces the gauge 
bosons to acquire mass in the SU(2) x [/(l)-symmetric electroweak theory. This method 
is an extension of the Abelian Higgs model we just studied. 

This mechanism postulates the existence of a doublet of complex scalar fields: 

$ = 4= ( t 1 1 2 ) (1-25) 

with real fields <fti, $2, <p3 and ^4. The most general S , t/(2)-symmetric renormalizable 
potential would be: 

V($) = y u 2 |$ t $| + A(|$ f $|) 2 (1.26) 
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with A > 0. 

Similar to the Abelian case, the vacuum acquires a v.e.v for /i 2 < 0; 



(<&> 



1 

71 



;i.27) 



shown in Figure (1.5). $ has the hypercharge Y$ = 1, and zero electromagnetic charge, 



Figure 1.5: The scalar potential with a v.e.v. leading to spontaneous symmetry breaking 

EDI. 



Q = I3 + y/2 = 0, which means that the electromagnetic U(l) symmetry is preserved. 
So, this mechanism yields the desired symmetry breaking SU(2)l x U(1)y — > U(1)em- 
The scalar field $ contributes to the Lagrangian with: 

L = (Dfty^Dpty - V($) (1.28) 

with the covariant derivative: 

D, = d, + i^-Bp + ■ (1.29) 
using the freedom to perform a gauge transformation on $, one could write 

*=4f ° 1 (1-30) 



V2\v + H 
We develop the theory around this minimum; 

6 + & 



$ = W VrrV-M (1.31) 

We substitute this $ in the Lagrangian, and study the second order terms to calculate 
the masses. The three fields are three Goldstone bosons which give mass to the three 
weak gauge fields. The spontaneous symmetry breaking also leads to a massive boson 
from the field H, the Higgs boson. 
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These mass terms are resulted form the following terms in the Lagrangian: 



WwiKf + (^m 2 ) 2 + (~9wW* + g B B,) 2 ] 



v 2 



Three massive vector bosons appear; 



The fourth vector field remain massless; 



with the weak mixing angle 9w] 



Wi = ^{W^iWl) , Z » = S^ + 9wK (L32) 
V2 V9w + 9b 



A , = 9w^ + 9bWI (L33) 
V9w + 9 2 b 



■ a 9w a 9b n OA , 

sin 6 W = — , cosg^ = — ; „ 9 (1-34) 

V + 9 B V9w + S'b 

The value for masses appears to be: 

M* / = ^g 2 w v 2 , Ml= l -{g 2 w + g 2 B )v 2 , M 7 = (1.35) 

Therefore, introducing the Brout-Englert-Higgs boson has managed to spontaneously 
break the electroweak symmetry and give mass to the gauge bosons while leaving the 
photon massless. 

Let's count the degrees of freedom in this case. We began with a complex scalar 
SU{2) L doublet $ with four degrees of freedom, a massless SU(2) L gauge field Wi with 
six degrees of freedom, and a massless U(l)y gauge field B with 2 degrees of freedom, 
making the total of 12. After the spontaneous symmetry breaking we have a physical real 
scalar field H (1 degree of freedom), massive W's and Z fields (9 degrees of freedom), and 
a massless photon (2 degrees of freedom), adding up to the same 12 degrees of freedom. 



1.3 Motivation for introducing more than one dou- 
blet 

There is no fundamental reason to assume that the scalar sector is minimal and contains 
only one Higgs doublet. The simplest extension of the Standard Model is the two-Higgs- 
doublet model (2HDM) [TT], with one extra Higgs doublet with the same quantum 
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numbers as the first one. There are many motivations for introducing two-Higgs doublet 
models, and they have been studied extensively. 

One of the motivations for introducing the second doublet is stemmed from the hi- 
erarchy of the Yukawa couplings in the third generation quarks, the top an the bottom 
quarks. In the Standard Model the masses of both these quarks comes from their cou- 
plings to the one Higgs doublet, which does not explain the ratio m t /mb ~ 35. However, 
this problem could be fixed if they received their masses from the two different doublets, 
providing the free parameters of the theory acquired appropriate values. 

One other reason to extend the scalar sector is to implement CP violation [TTJ in 
the model. In two-Higgs-doublet models the CP symmetry can be violated explicitly by 
choosing complex coefficients. 

Another motivation for introducing the second doublet is its application in super- 
symmetry [12]. The scalars appearing in supersymmetric theories are chiral multiplets, 
and their conjugates are of opposite chirality, therefore they cannot couple together. As 
a result, two Higgs doublets are needed to simultaneously give mass to the quarks with 
different charges, and therefore explaining the mass hierarchy in the fermion sector. 

Moreover, two-Higgs-doublet models are a possible source of Flavour Changing Neu- 
tral Currents (FCNC). These processes are strongly suppressed by experimental data, 
even though they do not seem to violate any fundamental law of nature. The Standard 
Model is compatible with experimental constraints on FCNC so far [T3J Q3] , with the re- 
markable exception of neutrino oscillations [15]. These processes are absent in the lepton 
sector, and in the quark sector they are prohibited. The simplest model that could ex- 
plain FCNC processes is the two-Higgs-doublet model. The Yukawa interactions arising 
from the second doublet lead to tree level FCNC's. 

Generating baryon asymmetry in the Universe of sufficient size is another motivation 
for introducing models with more than one doublet [16] , since they have a flexible scalar 
mass spectrum and they provide additional sources for CP violation. 

Models with more than one doublet, such as two-Higgs-doublet inert model |17j . 
could also accommodate dark matter candidates. In such models a Z 2 symmetry which 
remains unbroken, requires one of the doublets to have zero vacuum expectation value, 
and therefore not to couple to quarks and leptons. As a result this doublet, the inert 
Higgs doublet, could be a dark matter candidate. 



1.4 Two-Higgs-doublet model potential 

The Two-Higgs-Doublet Model (2HDM) is one of the simplest BSM extensions of the 
Higgs mechanism. In this model, one introduces two complex Higgs doublets 0i and (f>2 
with: 




(1.36) 
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The most general Higgs potential of 2HDM which is gauge invariant and renormalizable, 
is usually divided into quadratic and quartic parts: 



V 

v 2 
v 4 



v 2 + v 4 
1 



(1.37) 



+ 



Ai 

2 
1 



m ii(0i0i) + ^22(^2^2) + m 2 l2 {(j)\(j) 2 ) + m\ 2 "(0201/ 
^i0i) 2 + y (020 2 ) 2 + A 3 (0 t 1 0i)(0^ 2 ) + H^Mi&K 



+ 



Ae^) + XMfo) (0l0 2 ) + h.c. 



This potential has 14 free parameters [compare to only two real parameters in the poten- 
tial with one doublet (1.26)]; the real parameters m^, mf 2 , Ai, A 2 , A3, A 4 , and the complex 
parameters m\ 2 , A 5 , A 6 , A 7 . 



1.4.1 A simple case of 2HDM 



Here we briefly study a simple case of 2HDM, with Ai = A 2 and A 4 = A 5 



Af; — A7 



m 



12 



0. Therefore the potential simplifies to; 



V 



w u (010i) +m 22 (0 2 2 ) 



:i.38) 



+ ^i) 2 + ^(4fo) 2 + y (0100(402) 

In order to have a stable vacuum, the potential has to be bounded from below. It is 
enough to require the quartic part to be positive, therefore A > and A + A' > 0. We 
require the doublets to acquire an expectation value; 





Requiring the potential to have a minimum; 

9V - 





v 2 



1.39) 



90i 



dV 

d(p 2 



;i.4o) 



imposes the following conditions; 



^(-mfi + 2\v{ + \'v£) = 0, v 2 (-m 2 22 + 2\v\ + X'vf) = 



;i.4i) 



Therefore, we have four possibilities; 
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• The phase A with v\ — and v 2 = 

There is no symmetry breaking in this case, since the potential has one global 
minimum at the origin, with the conditions A > and A + A' > 0. To get the mass 
eigenstates, we develop the potential around the minimum; 

/ wf \ ( w+ \ 

fa= y h±m y fo= y J (1.42) 

which results in a mass spectrum of the form: 

2 2 

m\h a , Va ) = ^, ^ (1.43) 

with the extra condition m\ x ,m\ 2 < 0. 

• The phase B 1 with v± ^ and v 2 = 

In this case the potential has two global minima in 0i direction in the (0i,0 2 )- 
plane, at 0i = ±y. In this case only one of the Higgs fields, 0i, couples to the 
fermions. To study the potential, we excite the minimum; 

wt \ ( Wo 



(j>i = y ^i+fei+iqi J , <p2= y h 2 +i m J (1-44) 

Implementing this point in the potential, and vanishing the linear terms we get; 

v\ = ^ (1.45) 
which results in a mass spectrum of the form: 

m>?) = - 2Am ^ +A ' m " (1.46) 

2/r \ 2 -2Amo 2 + A'm?i . 

m \ha,T]a) = m u, tt (double degenerate) 

4A 

Positivity of the mass eigenstates requires; 

- ^22 + ^rn 2 u > (1.47) 

This condition is satisfied in three different cases; A > A' > 0, A' > A > and 
> A' > -A. 
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The phase B 2 with v± — and v 2 7^ 

This case is very similar to i?i case. Here the potential has two global minima in 
02 direction in the (0i, 02)-plane, at 02 — =t^j an d only one of the Higgs fields, 
02, couples to the fermions. Again to study the potential, we excite the minimum; 

wf \ ( W2 \ 

hi+im ) 02 — I ^2+^2+^2 (1-48) 

Implementing this point in the potential, and vanishing the linear terms we get; 

vl = (1.49) 

which results in a mass spectrum of the form: 

m>f) = -2Xml^X m > 2 (150) 

9/7 x 9 — 2Am?i + A'mno ni , N 

™> {ha,Va) = "^22? 77 (double degenerate) 

4A 

With the positivity of the mass eigenstates requiring; 

-m 2 u + jm 2 22 >0 (1.51) 

Similar to B x , here we have three cases depending on the values of A and A'. 

• The phase C with vi ^ and v 2 7^ 

The potential acquires four minima in this case, two in the 0i direction at 0i = 
±t>i/2, and two in the 02 direction at 02 = ±t> 2 /2. In this case the two Higgs fields 
can couple to fermions, with their different coefficients. We excite the minimum; 

_ / wf \ , _ ( wt \ 
0i — I , 02—1 v2+h 2 +im (1-52) 

Implementing this point in the potential, and vanishing the linear terms we get; 
9 4Am?i — 2A'mo 9 9 4Am 99 — 2A'm? 1 

4^-A° • ■* = w-a« (L53) 

which results in a 2 massless charged bosons, 2 massless neutral bosons and 2 
massive neutral bosons which are the eigenvalues of the following mass matrix: 



A \> D ~ m i2 i 3 \ 2 , 
4:^1^2 —4 r + 



(1.54) 
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Chapter 2 

N-Higgs-Doublet Models 



The electroweak symmetry breaking (EWSB) in the Standard Model is based on the 
Brout-Englert-Higgs mechanism (BEH). Many different variants of the BEH mechanism 
have been proposed so far, but it is not known which one is realized in the Nature yet. 
In the meantime, it is important to be aware of all essential possibilities that are realized 
in a chosen model. 

The non-minimal Higgs models have many free parameters, and it is important to 
analyse the model in its most general form, allowing for all possible degrees of freedom. 
This general analysis shows which phenomenological consequences are sensitive to the 
values of the parameters, which symmetries can arise in the model and how they are 
broken. After understanding the general structure of a model, one could restrict the 
model according to the existing experimental constraints. 

Unfortunately, such an exhaustive analysis is not possible in many non-minimal Higgs 
models. Even in the case of 2HDM [TTj [T5j [T9"] . the scalar potential cannot be minimized 
explicitly in the general case. Recently, a number of tools have been developed helping 
to understand the properties of the general 2HDM. These methods were based on the 
idea of the reparametrization symmetry, or basis-invariance, of the model: a unitary 
transformation between the Higgs doublets changes the parameters of the model, but 
leads to the same physical properties of the observable particles. 

This idea can be implemented via a tensorial formalism at the level of scalar fields 
[2"U| |2"TI [221 123] or through a geometric constructions in the space of gauge-invariant 
bilinears f2M [221 [251 [23 [2H] • In the latter case the formalism was extended to include 
non-unitary reparametrization transformations [291 EDI EI], which revealed interesting 
geometric properties of the 2HDM in the orbit space. 

It is a natural idea to extend these successful techniques to iV doublets. The general 
potential in multi-Higgs-doublet models is much richer than 2HDM, both at the level 
of scalar sector and Yukawa interactions (321 E31 El])- Some properties of the general 
NHDM were analysed in [261 EHl EH EH [39], with a special emphasis on CP- violation, 
[2T] HO]. However, a method to systematically explore all the possibilities offered with 
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iV doublets was still missing. 

Generalization from 2HDM to NHDM is straightforward in the tensorial formalism, 
however it is very difficult to translate tensorial invariants into physical observables. On 
the other hand, the geometric approach in the space of bilinears offers a more appealing 
treatment of the scalar potential, but the shape of the NHDM orbit space is rather 
complicated and has not been fully characterized so far. 

In this Chapter we study the algebraic and geometric properties of the NHDM orbit 
space. First we study the three distinct, but interrelated approaches to description of 



gauge orbits in the space of scalar fields in Section 2.1 Then we show a geometrical 



description of the orbit space in the general case of NHDM. In Section 2.2 we present 
the general scalar potential in NHDM, and end the Chapter with a detailed study of 
the orbit space of three-Higgs-doublet models, which is going to be used in the next 
Chapters. 



2.1 Orbit space of NHDM 

The scalar potential of the NHDM is constructed from gauge-invariant bilinear combi- 
nations of the Higgs fields. The space of these combinations can be described in three 
algebraically different but closely related ways: via representative Higgs doublets, the 
so-called i^-matrix, and a vector in the adjoint space. Here we describe and compare 
these three ways. 



2.1.1 Field space and gauge orbits 

The scalar part of the general NHDM has N complex Higgs doublets with electroweak 
isospin Y = 1/2: 

a=l,...,N (2.1) 

The dimension of the space of scalar fields is AN, with each doublet having 4 degrees of 
freedom (2 neutral and 2 charged). Since the scalar Lagrangian is electroweak-symmetric, 
we can perform any simultaneous intradoublet SU(2) x U(l) transformation inside all 
doublets without changing the Lagrangian. 

If we take a generic point in the Higgs space and apply all possible SU(2) x U(l) 
transformations, we will get a four-dimensional manifold called the (gauge) orbit. Thus, 
the entire 4iV-dimensional space of Higgs fields is naturally divided into non-intersecting 
orbits. The resulting set of orbits is a (4iV — 4)-dimensional manifold called the orbit 
space. 

When minimizing the Higgs potential, we look for vacuum expectation values of the 
Higgs fields (<p a ). Then, we can characterize each orbit by a specific representative point 
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in the Higgs space: 



\ , ( u 2 \ , ( u a e 



v i I ' ' ' \ v 2 e"^ I ' ' " V v a e l ^ a 



a>2 (2.2) 











0a = ( 





This point (and therefore, the entire orbit) is characterized by 4N — 4 real parameters: 
iV values of v a , N — 1 values of u a , a > 1, AT — 1 phases £ a , a > 1, and N — 2 phases r/ a , 
a > 2. 

A point with at least one u a ^ corresponds to the charge-breaking vacuum, in which 
the electroweak symmetry is broken completely and the photon acquires mass. For a 
neutral vacuum, we must set all u a = 0. Thus, the representative point of a generic 
neutral orbit is 

a>\ (2.3) 

This point is characterized by N parameters v a and N —1 phases £ a , making the dimen- 
sionality of the neutral orbit space equal to 2N — 1, which is 2N — 3 units less than the 
dimensionality of the entire orbit space. 

In addition to the electroweak SU(2) x U(l) transformations, which do not affect 
the parameters of the potential, one can consider the SU(N)h group of transformations 
that mix the doublets without affecting their intradoublet structure (index H stands 
for "horizontal" in contrast to the "vertical" transformations that mix the intradoublet 
structure). These transformations transform a given Higgs potential to another Higgs 
potential with different coefficients. Such a transformation is called a reparametrization 
transformation, (or a horizontal space transformation, or a Higgs-basis change). 

These transformations reparametrize the potential, but leave the physical observ- 
ables invariant [20, 21J. The same applies to anti-unitary transformations. Therefore, 
reparametrization transformations consist of all unitary and anti-unitary transformations 
U a b where 

<Pa^U ab b , <p a ~> Urifl ( 2 - 4 ) 



The antiunitary transformations are also known as generalized CP transformations 

32]. 

Reparametrization transformations link different gauge orbits, so if we pick up a 
generic point in the gauge orbit space and apply all possible reparametrization trans- 
formations, we reach many other points in the orbit space. This divides the orbit space 
into non-intersecting SU(N)-ovbits. 

Using reparametrization transformations, any neutral orbit parametrized by v a , £ a 
can be brought to a "canonical form" (also known as Higgs basis) 











0a = ( 





for a>l, v 2 = ^2v 2 a (2.5) 
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Depending on which doublet gets the non-zero value, there will be N equivalent canonical 
forms. 

Equivalently, any point in the charge-breaking orbit space can be brought to its own 
canonical form; 

There are N(N — l)/2 such canonical choices. To avoid double counting of equivalent 
canonical forms within such choices we can restrict v 2 > u 2 . The neutral orbit space 
corresponds to the limit u 2 — > 0. The other extremum, u 2 = v 2 , corresponds to a rather 
special space of "maximally charge-breaking" (MCB) vacua. 



2.1.2 K- matrix formalism 

One could define the electroweak-scalar bilinears as components of a complex hermitian 
N x N matrix [251 ISSl [27]: 

^ O 6 = (0l0a) (2.7) 

According to [261 the i^-matrix has the following properties; 

• It is a hermitian positive-semidefinite matrix. 

• Its rank is 2 for a charge-breaking vacuum and 1 for a neutral vacuum. This is 
because we are dealing with electroweak doublets and not higher multiplets. 

The fact that the .fT-matrix is a hermitian and rank 2 matrix, defines the dimensions 
of the neutral and charge breaking orbit space. It has at most two non-zero (and positive) 
eigenvalues, and the other N — 2 eigenvalues are zero. In the case of a neutral vacuum, 
the i^-matrix has only one non-zero (positive) eigenvalue and N — 1 zeros. 

If rank(f^) = 2, then among the N rows(columns) of K, there are at most two linearly 
independent rows(columns). Since this is an arbitrary choice, we require them to be the 
first and second rows (columns). All the other rows(columns) can be rewritten as linear 
combinations of lines 1 and 2. 

Since K is hermitian, these expansion coefficients are determined by the elements of 
the linearly independent lines [43J. Therefore, the number of algebraically independent 
gauge invariants K ab = 4>\4> a can be counted from 

K la = $ a (j) U a = l,...,N 

K Vt = <j>\<h, b = 2,...,N (2.8) 

where Ku and K22 are real while the rest are complex. Therefore the dimension of the 
charge breaking orbit space is 4iV — 4 [2iV from K\ a + 2(iV — 1) from K 2 b - 2 considering 
that Ku and K 2 2 are real]. 
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In the neutral orbit space rank(i^) = 1 and the first row/column is the linearly 
independent row/column, and therefore the dimension of the neutral orbit space is 2N — 1 
[2N from K\ a - 1 considering that K\\ is real]. 



2.1.3 Adjoint representation 

The adjoint representation of the group SU(N) H , is another way to look at the orbit 



space of the iV-Higgs-doublet model. Since the K-matrix (2.7) is a hermitian N x N 



matrix, it can be written as a linear combination of other hermitian N x N matrices; 



K = r ■ x j N ^_ 1} In + r 4 A, , % = 1, . . . , N 2 - 1 (2.9) 

where Aj are generators of SU(N) [43J. The coefficients tq and can be calculated from 
the .fT-matrix: 

ro = /Sr TrK ' r * =1 2 Tl [ KX ^ ( 2 - 10 ) 

The space of all possible real vectors r 1 ' = (r ,rj) is called the adjoint space. (The 
notation is similar to the Minkowski-space formalism developed for 2HDM in [271 129| 
EH ED].) 

The space of positive semidefinite i^-matrices with rank 2 is called the "i^-space". 



The mapping (2.9) from the K- space to the adjoint space r M defines a manifold called 
orbit space [37j . 

The orbit space has a complicated shape. In 2HDM the i^-matrix, a hermitian 2x2 
matrix, automatically satisfies ranki^ < 2, and the positive-semidefiniteness requires 

TrK > , (TriT) 2 - Tr[fT 2 ] > (2.11) 

which in the adjoint representation translates to 

r >0, r 2 >rf (2.12) 

which represents (the surface and the interior of) the forward light-cone in M 4 . 



With more than two doublets, there are extra conditions in addition to (2.11); 

n - 2 <^y<i (2.i3) 



2(JV-1) - \r„ 



The full conditions are formulated in 

The orbit space still lies inside the forward light-cone, but it occupies only a certain 
region inside it. 
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2.1.4 Geometric description of the orbit space 



Here we describe some geometric properties of the orbit space of NHDM in the adjoint 
space of vectors r M , which can be deduced from (2.13). 



From r 2 > rf, one could see that the orbit space has a conical shape. For convenience, 
we switch to the (N 2 — l)-dimensional space of normalized vectors rij = rj/r where the 
neutral orbit space lies on the surface of the unit sphere n 2 = 1, and the charge-breaking 
orbit space lies in a region inside it. 

The condition 

N-2 



2(N 



< nf < 1 



(2.14) 



shows that the orbit space is actually located inside a conical shell lying between two 
coaxial cones, the "forward light-cone" and a certain inner cone. Note that the orbit 
space does not occupy the whole space between these two cones. In the case of 2HDM, 
N = 2 the inner cone disappears, and the orbit space fills the entire cone. In NHDM, 



N > 2, this inner cone appears as shown in Figure 2.1 

















7, < 


7, 




(a) 


(b) 



Figure 2.1: The orbit space of 2HDM (a) versus the orbit space of NHDM (b) 



2.1.5 The orbit space in 3HDM 

In this Section we study the orbit space of the three-Higgs-doublet models, since certain 
examples of 3HDM will be presented in the next Chapters. 



From (2.9) we have 



K^ro—U + nXi (2.15) 
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The standard expression for the Gell-Mann matrices is the following: 



Ai 



A, 






(0 


—i 






)■ H 


i 








1. H 















f o 





>\ 




)• H 











|, *.= ( 




V -i 





0) 








1 











1 











-2 



A 7 = -% , A 8 = -= 1 (2.16) 



The explicit expressions for the field bilinears are: 

ro = i M - = m ± m ± m (2 . 17) 

V 3 v 3 

L [n , (4l0l) - (402) (0l0l) + (402) - 2(0^3) 

Tr[K Xi\ , r 3 - 



1 2 l *j - 2 , 2 ^ 

rj = Re(402) , r 2 = lm(0{0 2 ) , r 4 = Re(4^i) 

r 5 = Im(40i) , r 6 = Re(403) , r 7 = Im(4</>3) 

Note that A5 and therefore r 5 have opposite sign to the usual convention. This is done 
to expose the symmetry upon cyclic permutation. 

Recall that the orbit space of 3HDM sits between two coaxial shells 

\<nl<l , n l = ^ (2.18) 

As mentioned before, the orbit space has a complicated shape in the space of r's. How- 
ever, it is possible to project the orbit space onto a 3-dimensional space. We will see in 
Chapter [4] that the octahedral potential is linear in terms of the parameters x, y and z 
introduced as; 

x = n\ + n\ + nl , y = n\ + n\ + , z = n\ + n§ (2-19) 

It turns out that projecting the orbit space onto the 3-dimensional (x, y, z)-space, gives 
it a much simpler shape. 



The orbit space in the (x, y, z)-space 

The definitions of x, y and z satisfy 1/4 < x + y + z < I, meaning that the orbit space 
sits inside a 3-dimensional pyramid in the (x, y, z)-space. However it does not fill the 
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entire pyramid. We do not have the exact shape of the orbit space yet. Here we present 
the constraints that we have found on the borders of the orbit space inside the pyramid, 
which is represented in Figure (2.2). 




Figure 2.2: The pyramid in the (x, y, ,z)-space where the orbit space projection fits 
between the shown points. 

• x m i n and x max on X-axis (points F and A) 

On the XZ-plane y = and therefore 1/4 < x + z < 1, which represents a trape- 
zoid. With the extra condition z = on the X-axis the minimum and maximum 
values of x are respectively xp = 1/4 and xa = 1, which are realized through 

F , \ / wsin(7r/3) \ / usin(-7r/3) \ A .f°\f°\f° 



v J y v cos(7r/3) J y v cos(— 7r/3) J ' \ V J\ V J\ V 

• Zmin and z m ax on Z-axis (points G and B) 

On the XZ-plane 1/4 < x + z < 1, with the extra condition x = the minimum 
and maximum values of z on the Z-axis are respectively xq = 1/4 and xb — 1, 
which are realized through 

-C)(o)(o) ' -C)(o)(o 

• Vunique on 1^-axis (point -E) 

On the F-axis z = 0, which in terms of the doublets translates to \<f>\<j>i\ = v 2 . The 
most general case satisfying this condition is realized through v.e.v.'s of the form 

\ ( v sin a \ / v sin (3e iri3 
v cos ae im I y v cos /3e ir?4 
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With this v.e.v we calculate 

x — - [cos 2 a cos 2 r]2 + (cos a cos 8 cos(^4 — 772) + sin a sin /3 cos r] 3 ) 2 + cos 2 /3 cos 2 774] 

Since we are on the F-axis, we need x — 0. The following four sets of conditions 
vanish x: 

1: cosr/2 = 0, cos 774 = 0, sin a sin /3 cos 7/3 = — cos a cos 8 

2: coso; = 0, cos/3 = 0, cos 7/3 = 

3: cosa: = 0, cos 774 = 0, sin 8 cos 773 = 

4: cos77 2 = 0, cos/3 = 0, sinacos77 3 = 

Substituting any of these sets in y 

y — - [cos 2 a sin 2 772 + (cos a cos j3 sin(77 4 — 772) + sin a sin j3 sin 77 3 ) 2 + cos 2 (3 sin 2 774] 

results in y — 1/3, which proves that this point is the only point the orbit space 
touches the F-axis. This point is realized at point E with 



E 





iv 



v sina 

v cos a 



v sin (3e im 
v cos 8 



cos a cos 8 + sin a sin /? cos 773 = 



• Umax in YZ-plane (point C) 

In the FZ-plane a; = 0, which in terms of the doublets translates to Re((f)\(f)j) = 0. 
The most general case satisfying this condition is realized through v.e.v. 's of the 
form 

W W0 

^1 J V ^2 ; V 

With this v.e.v we calculate 



y = 



lm 2 (0l 



i'rjj 



3(7J 1 7J 2 ) i 



rl (vi 2 + tj 2 2 ) 2 

which has maximum value y max — 3/4 at V\ — v 2 . Therefore point C in the 
FZ-plane realizes a vertex of the orbit space with 



C 





iv 



• Umax in XY-plane (point D) 

In the XF-plane z — 0, which in terms of the doublets translates to |0|</>j| = v 2 . 
The most general neutral v.e.v. satisfying this condition is realized through 











ve 







ve 
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With this v.e.v we calculate 

y = - (sin 2 f] 2 + sin 2 (?73 - r] 2 ) + sin 2 1] 3 ) 

Trying to maximize y, one gets the condition sin(277 2 ) = — sin(277 3 ) = sin 2 (773 — 772) , 
which requires 772 = —773 or 772 = mr/2 ± 773, for 772 = nix/2 or 772 = nii/2 + 71/6. 
Substituting each of these four conditions in y, one gets y max = 3/4, which is 
realized in 

D: fo\f \ ( 



v J \ ve 2iw/3 ) \ ve- 2i7r/3 
The border line at y = 3/4 (edge CD) 

The border of the orbit space at y = 3/4 represents neutral vacua, therefore the 
upper components of the most general v.e.v. satisfying this condition are all zero. 
With y = 3/4, we also have the condition x + z = 1/4 which is satisfied in the 
most general case with 

\ / \ / 



V\ I \ v 2 e' l ' n J \ v 2 e %v 



With this v.e.v we calculate 

Im 2 (0l0,; 



with 



(l + 2r) 2 

2 



[2r sin 2 77 + r 2 sin 2 (27])] 



To maximize y, one requires cos(2r/) = — i. Substituting this value in y, results in 
y = 3/4 for any value of r, which means that x + z = 1/4 for any value of r. This 
condition represents the edge CD. 

The border line in X.Z-plane (edge AB) 

In the XZ-p\ane y = 0, which in terms of the doublets translates to lm(0|0j) = 0. 
The most general case satisfying this condition is realized through v.e.v. 's of the 
form 

'0\/0\/0 



, vi J \ v 2 J \ v 3 
With this v.e.v we calculate 

3{viv 2 ) 2 + 3{v 2 V3) 2 + (v 3 vi) 2 _ vf + v 2 + vl - v\v\ - v\v\ - v\v\ 



X 



(Vf + V 2 + V 2 ) 2 ' {v 2 +v 2 +v j)2 



For every value of Vi,v 2 and v 3 , x + z = 1, which represents the AB edge in Figure 
(pb. 
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The border line in YZ-pla.ne (edge BC) 

In the FZ-plane x — 0, the most general case satisfying this condition, He(<f)\(f)j] 
0, is realized through 

' o \ / \ / 



With this v.e.v we calculate 



Vi J \ iv 2 J V 



3(viv 2 ) _v\ + v\- vfv 2 

y (vi 2 + v 2 2 ) 2 ' Z Vi 2 + V 2 2 ) 2 

Introducing (j^j — r > an d rewriting y and z in terms of r 

3r 1 + r 2 — r 

y= (1 + r) 2 ' ^ = (1 + r) 2 

one could see that for any value of r, y(r) + ^(r) = 1, which represents the border 
line in the FZ-plane, the BC edge. 

The border line in X^K-plane (edge AD) 

In the Xy-plane z = 0, the most general v.e.v. satisfying this condition, \4>\4>i\ = 
v 2 , is realized through 

' \ / \ / 



v I \ ve"' IA I \ ve %m 



With this v.e.v we calculate 



y = - (sin 2 772 + sin 2 (?73 - rj 2 ) + sin 2 r/ 3 ) , x = - (cos 2 r] 2 + cos 2 (r] 3 - r) 2 ) + cos 2 773) 
For every value of rj 2 and 773, x + y — 1, which represents the AD edge in Figure 



(2.3). 



So it seems that the orbit space lies inside the shape shown in Figure (2.3). 



2.2 The scalar potential in NHDM 

In this Section we show how the scalar potential of multi-Higgs-doublet models is pre- 
sented in the orbit space^j The generic Higgs potential in NHDM in terms of the doublets 
is written in a tensorial form 1201 121] : 



V = YMlfo) + Z, 65d (0t0 fe )(0t 0d) (2.20) 



lr rhe idea to switch to the orbit space in order to simplify the task of a group-invariant potential 
minimization is rather old, see 44 and references therein. 
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Figure 2.3: The borders of the orbit space of 3HDM in the (x, y, z)-space. 

where all indices run from 1 to N. The potential is constructed from iV 2 bilinears 
{<f>a<j>b)- Coefficients in the quadratic and quartic parts of the potential are grouped into 
components of tensors Y & b and Z ahsd , respectively; there are iV 2 independent components 
in Y and N 2 (N 2 + l)/2 independent components in Z. 

In the adjoint space the Higgs potential has the following form: 



where the tensor Y n b is written as the scalar Mo and the vector Mi in the adjoint space. 
The number of free parameters in Mo and the vector Mj is [1 + (iV 2 — 1) =]N 2 . The tensor 
Zabcd is represented as the scalar A o, vector A j and symmetric tensor Ay in the adjoint 
space. The number of free parameters here is [l+(N 2 -l)+N 2 (N 2 -l)/2 =]N 2 (N 2 +l)/2. 

Using the adjoint space r M makes the analysis of the Higgs potential in NHDM much 
simpler. It is easy to formulate and prove conditions for existence of a symmetry. We 
will present examples of this in Chapter [3] and [4} 




) (2-21) 
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Chapter 3 

Frustrated symmetries in NHDM 



In this Chapter we introduce a peculiar feature that arises in models with more than two 
Higgs doublets. We show that there exists certain explicit symmetries of the potential, 
which are necessarily broken whenever the electroweak symmetry breaking takes place. 
We term them "frustrated symmetries" because of their similarity to the phenomenon 
of geometric frustration in condensed-matter physics. We think that such symmetries 
can be an interesting option when building models with desired properties beyond the 
Standard model. 



The structure of the Chapter is the following; In Section 3A_ we explain how frustrated 
symmetries arise in NHDM with N > 2. In Section 3.2 we give several examples of such 
symmetries in 3HDM. The results of this Chapter are published in 



3.1 General observations 

Suppose we have constructed a Higgs potential with some explicit symmetry. That is, 
we have found coefficients Mj, A 0i and Ay in the potential: 

V = -M r - Min + ^A 00 ro + A 0i r 7i + ^A^nr, (3.1) 

such that the potential is left invariant under group G of transformation between dou- 
blets <f> a . We will restrict ourselves to reparametrization transformations which are either 
unitary or antiunitary transformations, and leave the Higgs kinetic term invariant. This 
means that vq is invariant under G, while the components are transformed by orthog- 
onal transformations. Note that: 




where Aj's are generators of SU(N). 
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In particular, it is possible to devise such a symmetry group G of transformations of 

ij 



(C) 

doublets, which can be implemented via a non-trivial tensor A) ? - , but not in a vector 
Mj or A Qi . Let's study a potential with such symmetry; 



= -M r + -A 00 r 2 + -Afnr, (3.3) 



2 2 

After EWSB, the doublets acquire some vacuum expectation values; 



1 



<« = 7ii":j (3 ' 4) 

Therefore (ro) proportional to (0a0 a ) has a non-zero value and (r^) is proportional 

If all (rj) = 0, then for any pair of doublets (4>\<fib) = 0. Therefore the potential, 
depending only on r 7^ which is invariant under G, conserves the symmetry G after 
EWSB, and the vacuum remains invariant under G. 

If at least one component (r») 7^ 0, then the vacuum spontaneously breaks the sym- 
metry group G, either completely or down to a proper subgroup, simply because by 
construction no non-trivial vector (r^) can be G-symmetric. 

Thus, in order to see if a given symmetry has a chance to be conserved in a non-trivial 
vacuum, we need to check whether or not the point n = (r) / (r ) = belongs to the 
orbit space of the model. 

In Chapter [2] we showed that the point n = belongs to the orbit space only in the 
two-Higgs-doublet model. Therefore, in 2HDM any symmetry imposed on the potential 
can, in principle, be conserved in the vacuum state provided we have chosen appropriate 
coefficients of the potential. Indeed, if one acts on the state; 

■ <*> = 7i(S) (3 ' 5) 

with any (anti)unitary transformation which mixes the doublets, one arrives at the same 
state up to an electroweak transformation. That is, the corresponding gauge orbit is 
G-symmetric. 

In NHDM with N > 2, the point n = cannot be realized through doublets, since 

A ~ ~ <n 2 <l (3.6) 



to 



2(N- 1) 



Therefore, the minimum of the G-symmetric potential (3.3) unavoidably breaks the 
explicit symmetry after EWSB. 

We call such a symmetry a frustrated symmetry for the following reason: setting 
all rii = is equivalent to setting all products {(jr a (pb) to zero while keeping all the norms 



a'ra 



non-zero and equal. This is impossible to achieve with more than two doublets, 
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simply because we have too little freedom of where to place vacuum expectation values 
inside all three doublets. 



Even if the first two doublets are chosen as in (3.5), the third doublet of the same 
norm will have a non-zero product with the first or the second doublet. 



<*> = ;ji(°) • <« = 7i(S) ' w = ^(0 (37) 

In other words, although we can "optimize" the v.e.v.'s in any pair of doublets (equal 
norms, zero product), these optimal configurations are mutually incompatible in three or 
more doublets and cannot be satisfied simultaneously. This is precisely what frustration 
in condensed- matter physics is about. 



3.2 Examples of frustrated symmetries in 3HDM 

In this section we present examples of frustrated symmetries we have found in the three- 
Higgs-doublet model. 



3.2.1 First example 

The simplest case is given by the potential which is symmetric under any £77(3) rotation 
and mixes the three doublets. It corresponds to the exceptional case A^- = 0: 

V = -M r + A r 2 Q 

v = -^(4^ + 4^ + 4^ + ^(4^ + 4^ + 4^ (3.8) 

Clearly, no non-trivial vacuum can be symmetric under the entire SU(3) group, even up 
to EW transformations. We choose: 

(&> = 4 ( H ) (3-9) 



Let's study the excitation around this point: 



-j^(v + h a + irja) 
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And rewrite the potential in terms of the new parameters: 



V 



Mo , n N , A 0/ _ 3 



+ 
+ 
+ 
+ 



2VI (2 " ) + Ii (12 " 

Mn An , 

+ -^-(10v 2 ) 



2v/3 12 

M A 2 

+ — (6u J 



(/ii + h 2 + h 3 ) 
(h\ + h 2 2 + h 2 3 ) 

(v! + vl + vl) 



2^3 12 

[h\h 2 + h 2 h 3 + h 3 hi) 
M , A 



> 2 



12v z 



V3 12 
+ higher order terms 



(wfw 1 + W2W2 + w^w 3 



(3.11) 



From vanishing linear terms, one gets v 2 = ^j^, which results in the following mass 
spectrum; 



m 2 (r] a ) 
m 2 (w±) 

m 2 (h a ) 








M 
3^ 



double degenerate 



(3.12) 



The massless bosons are of course expected, since the original symmetry of the potential 
is continuous (Goldstone Theorem). 



3.2.2 Second example 

The second example is given by the " tetrahedral" symmetry in 3HDM, defined by the 
following scalar potential: 



V = -M r + Aor 2 + Ai(r 2 + rl + rl)+A 2 (r 2 2 +r 2 + r 2 ) 
+A 3 (r 2 + r 2 s ) + A 4 (r!r 2 - r 4 r 5 + r 6 r 7 ) 



(3.13) 
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In terms of doublets, this potential has the form: 



V 



Mi 



Ao 
3 

Ai 
3 

+A 2 



[fa + fa 2 fa + fa 3 fa) 



ifa) + (4<p 2 ) + (4fa) 



3^1 ) 



+A 3 
+A 4 



(^> + i0i) 2 + (4fa) 2 + (4fa) 2 - {4faMfa) - (4faMfa) - (4faMfa) 

(Re0 t 1 2 ) 2 + (Re0^ 3 ) 2 + (Re^ t - i > 2 



(3.14) 



(Im<M0 2 ) 2 + (lm4fa) 2 + (Im^fa) 
(Re0 t 1 2 )(lm0 t 1 2 ) + (Re4fa)(lm4fa) + (Re4fa)(Im4 3 fa) 



This potential is symmetric under the chiral tetrahedral group (Appendix B). 



Note that the absence of any term linear in r, in (3.13) is due to the fact that no 
such term respects the tetrahedral symmetry. 

For Ao < we have a vacuum of the form: 



(fa) 





V2 V v 



(3.15) 



And study the excitations around this point: 



w~ 



^(v + h a + irja) 
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(3.16) 



The potential has the following form in terms of the new parameters: 



V 



(/*! + h 2 + h 3 ) 



+ 



" Mo + ^(10^ 2 ) + ^(4t; 2 ) + ^(2^) 



2x/3 12 



12 



(hi + hj + h 2 3 ) 



B 



+ 



— 1=- H 6f H 2t> 

2v/3 12 V ; 4 V ; 



a 



+ 



- Mo + ^(12^ 



y/3 12 



(Wi + Wi + w 2 + w 2 +w 3 + w 3 ) 



D 



+ 



(Wi w 2 + + w 2 W 3 + + W 3 Wi 



H 



+ 



(wfw 2 + w^w 3 + w^w 1 ) 



2E 



+ 



(/ll/l 2 + /i2^3 + h^hx) 



2F 



+ 



(?7i?72 + ^3 + mm) 



2G 



+ 



A 



> 2 ) 



(^-l?72 - Ml + ^3 - ^3^72 + Ml - hffy) (3.17) 
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From the vanishing linear terms v 2 



A 



,A + Ai 



Mp 

v^Ao+Aa)' 

B 



which simplifies the coefficients in (3.17) 



-A 2 + A 3 2 
)v 



C 



-A 2 )v 2 



2A 2 + iA 4 2 

H = ; )V 



d = ( 2A *- iA *y 

2E = ( 2A °- Al + 3A2 )^ 



2F 



.A 2 -A 3 2 



2G 



4n 2 



2 y ' v 2 

These coefficients result in a mass spectrum of the form; 

2/ ±x 6|A 2 | ± v / 3A 4 2 
m (w a J : ; v 



4 



(3.18) 



m 2 (h a ) 
m 2 (r] a ) 



v 

T 



2A 1 + |A 2 |+3A 3 ±VA 

2A 1 + |A 2 |+3A 3 ±v / A 
with A = 4A 2 + 2|A 2 | 2 + 9A 2 - 20Ai|A 2 | - 12AxA 3 + 30|A 2 |A 3 + 3A 

The positivity of mass eigenstates imposes constraints on A's; 

-6|A 2 | < V3A 4 < 6|A 2 | 
|A 2 | 2 + 3A 2 < 24(AxA 3 - |A 2 |A 3 + A 1 |A 2 |) 



(3.19) 



3.2.3 Third example 

The third example of frustrated symmetry is the "octahedral" symmetry in 3HDM, 
whose potential is given by; 

V = -M r + A r 2 + A x (r 2 + r\ + r 2 ) + A 2 (r 2 2 + r 2 + r 2 ) + A 3 (r 2 + r 8 2 ) (3.20) 

This potential is symmetric under the octahedral group (Appendix B). 

If A 2 < 0, the following v.e.v.'s realize the global minima of this potential; 



As usual we check the excitations of the fields around this point: 

(0a) = 



(3.21) 



-j^(v + h a + irja) 
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(3.22) 



then rewrite the potential in terms of the new fields: 



V 



2x/3 12 V ; 4 V ; 

A 





+ h 2 + h 3 ) 



+ 



" M ° + ^(10v 2 ) + ^(2v 2 ) + ^(Av 2 ) 



2x/3 12 



iWn An / on A3 . ox 

+ — y= + — (6v 2 ) + — (2v 2 ) 



+ 




12 



(hf + h 2 2 + hi) 



+ -^ + ^l2v 2 ) 



(w^w 1 + w 2 w 2 + w 3 w 3 



w 1 w 2 + w 2 w 3 + w 3 + W{W 2 + w 2 w 3 + w 3 w 1 ) 



+ 



(/ii/i 2 + ^2^3 + h 3 hi) 



+ 



2F 



(3.23) 



With — ^3(X^+X^)' we s i m P^fy the coefficients 



,A + A 2 2 
v o > 



2E 



C = (-A 2 )v 2 

-y 



,2A + 3A 2 -A U 2 



,-A 2 + A 3 , o 

5 = 2 —^) v 2 

y 2 ' 

,A 2 - A 3 , o 
2F = (— -)v 2 
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The eigenvalues of the mass matrix are: 

3|A S 



m (w a ) : 

m 2 (h a ) : 
m 2 (rj a ) : 



l 2 2 



(double degenerate) 



|A 2 



Ai 2 

V 



6A -8|A 2 | -2Ax 2 

V 

3 3 
(3|A 2 | + 3A 3 )t> 2 (double degenerate) 



(3.24) 



Positivity of the eigenvalues of the mass matrix puts the following constraints on A's: 



A 2 < , Ai + |A 2 | > , A - |A 2 | > 
A 3 >0 , A 3 + |A 2 |>0 



(3.25) 



We note a remarkable phenomenological feature of this model: it is 2HDM-like. 
Due to remaining symmetry, it exhibits certain degeneracy in the mass spectrum of the 
physical Higgs bosons, yielding just one mass for both charged Higgs bosons and three 
different masses for the neutral ones, which precisely mimics the typical Higgs spectrum 
of 2HDM. This model is studied in detail in Chapter |4} 



3.2.4 Forth example 

The forth example is of an Abelian frustrated symmetry in 3HDM; a Z% x Z% symmetric 
potential, which is given by; 



V 



+A 3 



(400 + (4^) + (4<j> 3 ) 



Ao 
3 

+ A 2 



^ 2 | 2 + l40 3 | 2 



{4Mmh) + (020 3 ) (0200 + (4^i)(m 



+ h.c. 



r 



(3.26) 



We again pick a vacuum point: 





V2 V v 



(3.27) 



And check the excitations of the fields around this point: 



}=(v + h a + irja) 



V2 



(3.28) 
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Then rewrite the potential in terms of the new fields: 



V 



-M 
x/3 



» + (A + Ai + A 2 + A 3 + A 3 *)(v 3 ) 



(hi + h 2 + h 3 ) 



+ 



A 



-M A 2 A x 2 A 2 . 2 A 3 + A 3 * 2 ^ 
+ To ( ) + ^( 6u ) + ^( 2w ) + ; ( u ) 



2^3 12 



(hi + hj + hi) 



+ 



-M A 2 A x 2 A 2 . 2 A 3 + A 3 * 2 
: (6v ) + — (2v ) + — (2v ) (v ) 



2x/3 12 



(r,j + V 2 2 + rjj) 



+ 



c 



-Mi 



+ A v 2 + A lV 2 



+ 



D 



A 2,^-3 + A 3 * 2 

A 2 v H [v ) 



(WiW 1 + W 2 W 2 + w 3 w 3 



(w 1 W2 + w 2 w 3 + w 3 + wfw 2 + W2 w 3 + W3" ) 



+ 



2E 



— (8v ) + A 2 v + (5v ) 



(hih 2 + h 2 h 3 + h 3 h\ 



+ 



2F 



A 3 + A 3 * 2 ^ 
a ( v ) 



+ 



G 



■ A3 - A3* 2 

I (v z ) 



(mm + mm + mm) 



(him + h 2 m + h 3 m - him - h 2m - h 3 m) 



(3.29) 
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With v 2 = + Ai +l 2 + Aa + Aa * we simplify the coefficients of the potential: 



2A + 2A i -A3-A 3 \ 2 

A = ( o ) V 



C = (-A 2 - A 3 - A 3 > 2 

_ .8A + 12A 2 + 15A 3 + 15A 3 * 2 
2£/ = — )v 



12 



G = i{ 



A 3 -A 3 \ 2 



D = { 
2F = 



-3A 3 -3A 3 \ 2 



2A 2 + A 3 + A 3 * 2 
)u 

), 2 



2 

A 3 + A 3 * 



The eigenvalues of the mass matrix are: 



2 -4A 2 -3A 3 -3A 3 2 

A 2 t> , -y 



(3.30) 



m 2 (h a ) 



m 2 (r] a ) 



with 



12 
18 
12 



(16A + 24Ai - 24A 2 - 27A 3 - 27A 3 *) 
(20A + 12Ai + 12A 2 + 9A 3 + 9A 3 *) 



3A + 12Ai - 6A 2 - 21A 3 - 21A 3 * + VA 



v 
12 



7v_ 2 
' 4 

2 r 



(A 3 + A 3 *), - W 2 (A 3 + A 3 *) 



8A + 12Ai - 6A 2 - 21A 3 - 21A 3 * - 



A = (8A + 12Ai - 12A 2 - 3A 3 - 3A 3 *) 2 + 



9(3^3-15) 



(A 3 -A 3 * 
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Chapter 4 

3HDM with octahedral symmetry 



In this Chapter we focus on a specific discrete symmetry group that can be imposed 
on the scalar potential of the three-Higgs-doublet model: the full (achiral) octahedral 
symmetry group Oh- We see two main goals in this study. 

First; Oh is among the largest realizable finite symmetry groups which can be imposed 
on the scalar potential in 3HDM. The word "realizable" stresses that when imposing such 
a symmetry group, we obtain a potential that is symmetric exactly under this and not 
any larger group. It is therefore interesting to check what are the phenomenological 
consequences of such a high symmetry of the potential. 

Second; the octahedral 3HDM serves as a good illustration of the power of geometric 
and group-theoretic methods which were developed in [43J for the space of Higgs bilinears 
and which are further developed in this Chapter. We believe that these methods can be 
used in the analysis of other similar models. 



The structure of the Chapter is as follows: In Section 4.1 we introduce the octahedral 



3HDM potential and describe its symmetry group and symmetry breaking possibilities. 



In Section 4.2 we show the main geometric idea in order to minimize the potential, with 



a simplified example. Then we apply the method to the full 3-dimensional orbit space 



and derive positivity conditions. In Section 4.3 we perform a phenomenological analysis 



on each possible type of symmetry breaking, and in Appendix [C] we study a specific 
minimum in more detail. 



4.1 The scalar potential 

The Octahedral 3HDM potential is defined by 

V = -M r + A r 2 + A 1 (rl + rj + r 2 6 ) + A 2 (r| + r\ + r 2 ) + A 3 (r 2 3 + r|) (4.1) 
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To write the potential in terms of the doublets, we remind the reader of the explicit 
expressions for the field bilinears in 3HDM: 



r 

n 

r 5 



+ 



+ 



Re{<f>\<t> 2 ) 



r-2 
r 6 



_ (4m + - 2(0^3) 
rs " 27! 

Im(0{0 2 ) , r 4 = Re(030i) 
Re(0 2 03) , r 7 = lm(0 2 3 ) 



Therefore the potential has the following form: 

Ao 

+Ai 



V 



+ 



101 + 0202 + 0303 



+A 2 
A 3 



(4.2) 



(4.3) 



(Re^) 2 + (Re0 2 3 ) 2 + (Re0 3 0i) 

(4^1? + (4 2 <p2? + (0 3 3 ) 2 " (0l0l)(0502) " (0502)(0303) - (0 3 3 )(0 t 1 01 



This potential is symmetric under full achiral group Oh (Appendix B) 



4.1.1 Symmetry breaking patterns 

It seems that in 3HDM a sufficiently symmetric group G, such as the octahedral sym- 
metry, can not break down completely after EWSB, since it will create so many minima 
[4*6] . So, it breaks down to certain subgroups of G, which narrows down the list of pos- 
sible patterns of v.e.v.s. Several such patterns can be simply guessed. Whether they are 
actually the minima of the potential for some values of the parameters, will be studied 
later; with this group-theoretic argument we can only list patterns which have a chance 
to be minima. 

Below we list several such patterns. We will focus on neutral minima only, and 
therefore we present only the lower components of the doublets. 

• (v, v, v) . This minimum is invariant under permutations of the doublets and under 
CP; the symmetry group is D 3 x Z 2 of order 12. Applying transformations which 
do not leave it invariant, we obtain four minima of this type: (v,v,v), (v,v, —v), 
(v, —v, v), and (v, —v, —v). 

• (v, 0, 0). This point is invariant under sign flips, under CP and under transposi- 
tion of the last two doublets, which form the symmetry group -D 4 x Z 2 of order 16. 
The three minima of this type are (v, 0, 0), (0, v, 0), (0, 0, v). 
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• (v, v, 0). This point is symmetric under the sing flip of the third doublet, permu- 
tation of the first and second doublet, CP transformation and their combinations. 
The symmetry group is {Z 2 ) 3 of order 8. The six minima of this type are (v, ±t>, 0), 
(v,0,±v), (P,v,±v). 

• (v, iv, 0). This point is invariant under the sign flip of the third doublet, the CP 
transformation compensated by the sign flip of the second doublet, exchange of 
the first two doublets compensated by sign flip of the second doublet, and their 
products. The symmetry group is Dh 4 of order 8. The six minima are (v, ±iv, 0), 
(v,0,±iv), (0,v,±iv). 

• (d, ^e 1 ^ 3 , De""'' 3 ) and (v, ue 2Mr / 3 , t;e _2 " r / 3 ). The symmetry group is S s and 
is made of permutations of the three doublets compensated by appropriate sign 
flips and/or CP. There are six distinct minima of the first type and two distinct 
minima of the second type. 



This seems to exhaust all sufficiently symmetric points. In the next Section we will 
find which ones can actually be minima of the potential. 



4.2 Minimization of the scalar potential 



The octahedral 3HDM potential (4.1) looks deceptively simple. If one plunges into 
straightforward minimization of the potential, it becomes a challenging task to recon- 
struct the full phase diagram of the model: that is, to find the positions (and even types) 
of the minima for arbitrary values of Aj. Here we develop an efficient geometric method 
to solve this problem. Before we proceed with the solution, let us illustrate this method 
with a simple toy model. 



4.2.1 Toy model 

Suppose only real- valued fluctuations in all three doublets are allowed (r 2 = r 5 = r 7 = 0). 
Then the potential simplifies to 



V = -M r + r 2 (A + A lX + A 3 z 
with 



2 , 2 , 2 

x = n x + n 4 + n 6 , 



2 , 2 

z = n 3 + n g 



(4.4) 



Recall from Section |2.1.5| that in 3HDM 1/4 < nf < 1, therefore 1/4 < x + z < 1, which 
defines a trapezoid shown in Figure (4.1). 
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% 1 x 

Figure 4.1: The orbit space of the toy model projected on the (x, z)-plane. 



The four vertices of this trapezoid correspond to the following patterns of the three 
doublets: 



x — , z — 1 : 
x = 1 , z = : 
x = 0, 2 = 1/4: 
x = 1/4, 2 = : 




The first two points correspond to the neutral vacuum, and they were already listed in 
the previous section, while the last two are charge-breaking minima. One can also show 
that all points inside this trapezoid are realizable by fields (in other words, when the full 
orbit space is projected on the (x, z)-plane, it covers the entire trapezoid). 

Positivity conditions 

The positivity conditions are restrictions on the parameters which guarantee that the 
potential is bounded from below for all values of the Higgs fields. It is known that it is 
necessary and sufficient to require that the quartic part of the potential is positive inside 
the orbit space. 
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Now comes the crucial point. Since the potential (4.4) is a linear function in x and z 



defined inside this trapezoid and therefore convex, it is sufficient to write the positivity 



conditions at the four vertices (4.5); then they will be automatically satisfied at all points 



inside the trapezoid. Thus, we arrive at the following positivity conditions for our toy 
model: 

A + Ai > , A + A 3 > , 4A + A x > , 4A + A 3 > (4.6) 

In general, none of these inequalities can be removed because A can be positive or 
negative. 

Now let us show how the global minima can be found for different values of the 
parameters. Since the potential is a linear function of x and z, we can introduce on the 
(x, z)-plane the direction of steepest descent (DSD), which is given by vector (— Ai, —A3). 
Then, we need to find what points of the trapezoid lie farthest along this direction; such 
points give the (x, z) coordinates of the global minima, that is, patterns of v.e.v.s of the 
doublets. The last step is to find the overall v.e.v. normalization v. 

By checking all possible directions of steepest descent, one finds situations where 



vertices or edges of the trapezoid correspond to the global minimum (see Figure(4.2)). 




edge AB 




► Ai 



Figure 4.2: (a) Various possible directions of steepest descent and the corresponding 
minima, (b) Phase diagram of the corresponding potential. 



If DSD is parallel to direction 1, which happens at Ai = A3 > 0, then the entire 
edge AB corresponds to the global minimum. If DSD is between directions 1 and 2 
(Ai > A 3 > 0), point A is the minimum. When DSD is along direction 2 (Ai > 
0, A 3 = 0), all the points on the AD edge are minima. When DSD rotates further from 
direction 2 to direction 4, point D becomes the minimum. When DSD reaches direction 
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4 (Ai = A3 < 0), then all the points on the DC edge become minima of the potential. 
And so on. In this way we can construct the full phase diagram of the model in the Aj 
parameter space. 



4.3 Different minima 



In this section we analyse the potential for each of the minima found in Section AJ2_ We 
focus only on the four possible points corresponding to the neutral vacuum. 



4.3.1 The minimum (v,v,v) 

In order for this point to be a minimum, it is necessary to have Ai < 0. The other 
parameters, A 2 ,A 3 , can be of any sign, but if negative, they cannot be lower than a 
certain limit. Explicit expansion of the potential at this point gives the following scalar 
mass spectrum: 

3 

m2 (ha) '■ T^^" 1 ^ 2 (double degenerate) 

m2 (^a) : (|Ai|+A 3 )f 2 (double degenerate) 
3 

-(| Ai| + A 2 )t> 2 (double degenerate) 

2(A -|A 1 |)t; 2 (4.7) 

The positivity of the mass eigenstates gives extra conditions on A's: 

Ai<0, A 3 + |Ax| > 0, Ao - |Ai| > 0, A 2 + |A!|>0 (4.8) 

The remarkable feature of this spectrum is that it is 2HDM-like, namely, it has only 
one massive charged scalar and three massive neutral scalars. We study this interesting 
case in more detail in Appendix C. 



4.3.2 The minimum (v, 0, 0) 



In order for this point to be the minimum, A3 must be negative. Checking the excitations 

V3M 
A -|A 3 |' 



around the minima, will result is v 2 = A . , and the mass spectrum is: 



m 2 (h^) : "^T^ 2 (double degenerate) 



2(A - |A 3 J^ 2 



3 

Ai + |A 3! 2 



A 2 + I A, 



v (double degenerate) 
-v 2 (double degenerate) (4.9) 
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With this breaking pattern, the spectrum is again 2HDM-like. Positivity of the mass 
eigenstates puts extra condition on A's: 



A > 0, A - |A 3 | > 0, Ax + |A 3 | > 0, A 2 + |A 3 | > 



(4.10) 



4.3.3 The minimum (v,iv,0) 



6 M 



This point appears as a minimum when Ai > 0,A2,A3 < 0. For v = _ 3 | A 
the mass spectrum has the following form: 



m\ht) : |A 2 |^ 



|A 2 | + |A 



3 2 
V 



m\hl) : 



Ai + |A 3 | 



v 2 (double degenerate) 



3Ai + 3|A 2 | 2 2 4A -3|A 2 | - |A 3 | 2 
v, (|A 2 | - \A 3 \)v , v 



(4.11) 



There appear extra constraints on the A's from the positivity of the mass eigenstates 
here: 



|A 2 | + |A 3 |>0, Ai + |A 3 | > 0, Ai + |A 2 | > 
|A 2 |-|A 3 |>0, 4A -3|A 2 | - |A 3 | > 



(4.12) 



4.3.4 The minimum (=bu, ve t7r / 3 , ve i7r / 3 ) 

For this point to be minimum, one requires A 2 < 0, Ai < |A 2 | and Ai < A 3 . The mass 
spectrum has the following form 



m 2 {h° a ) : (4A -3|A 2 | + A 1 )^ 
"3(|A 2 |+Ai) 



3|A 2 |^ 2 3(jA 2 j-A 1 ) u 2 



+ A 3 - Ai ± 



v 2 (double degenerate) 



with A = 



3(|A 2 | +A0 



+ (A 3 - A,f 



(4.13) 



This spectrum shows a bizarre pattern: the charged Higgs bosons are non- degenerate, 
while the neutrals are 2HDM-like. 
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Chapter 5 

Abelian symmetries in NHDM 



5.1 Symmetries in NHDM 

One of the interesting branches of study in non-minimal Higgs sectors is the possible 
symmetries one can impose on the potential. In particular, three topics are of importance: 

1. Which symmetries can a general NHDM potential have? 

2. What are the possible ways to spontaneously break these symmetries? 

3. How can these symmetries be extended into the fermion sector? 

In this Chapter we explore the first question. When studying the symmetries of 
the potential, one focuses on the reparametrization transformations, which mix different 
doublets but leave the intradoublet structure untouched. These transformations keep 
the kinetic term invariant, and are either unitary or antiunitary. 

In this Chapter we focus on the unitary transformations. Such transformations form 
the groups U(N). We are looking for the physically distinct transformations, therefore 
we leave out the common phase rotations since they act trivially on the potential. The 
overall phase factor multiplication which is taken into account by the gauge group U(l)y 
should be left out. This leaves us with the special unitary group SU(N) which has to 
be factor-grouped by its center Z N . We call the resulting group PSU(N). We construct 
this PSU(N) for every NHDM, and look for the maximal Abelian subgroups of it, and 
then explore all its realizable subgroups. 

There are several maximal Abelian subgroups in PSU(N), but they are all conjugate 
to each other. Therefore we could pick one maximal Abelian subgroup and exhaust its 



subgroups. Figure (5.1) shows the tree of subgroups in NHDM. 

When searching for the symmetry group which can be implemented in the scalar 
sector of a non-minimal Higgs model, one must distinguish between realizable and non- 
realizable groups, [471 . If it is possible to write a potential which is symmetric under a 
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group G but not symmetric under a larger symmetry group containing G, we call G a 
realizable group. A non-realizable group, on the contrary, automatically leads to a larger 
symmetry group of the potential. The true symmetry properties of the potentials are 
reflected in realizable groups. 

For example, a potential which depends on the first doublet only via (<f>\<j>i) is ob- 
viously symmetric under the cyclic group Z„ of discrete phase rotations of this doublet 
for any n. However, these Z n 's have no interest on their own because they trivially arise 
as subgroups of the larger symmetry group of this potential U(l) describing arbitrary 
phase rotations. It is this U(l), not its individual subgroups Z n , which has a chance to 
be the realizable symmetry group of this potential. 

Therefore, if one focuses on the scalar sector of the model and aims at classifying the 
possible symmetries which can be implemented in the scalar potential, it is natural to 
restrict one's choice to the realizable symmetry groups. 

In this Chapter we introduce the group of reparametrization transformations for 
NHDM. Next we develop the strategy to find all unitary Abelian realizable symmetry 
groups in NHDM. Using this strategy we present the full list of unitary Abelian realizable 
symmetries in 3HDM and 4HDM. In the attempt to try to derive the full list of unitary 
Abelian realizable symmetries in NHDM we prove several statements for a general NHDM 
potential. The results of this Chapter is published in |J8]. 
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5.2 Scalar sector of the iV-Higgs doublet model 



5.2.1 The group of reparametrization transformations 

When discussing symmetries of the potential, we focus on the reparametrization trans- 
formations, which are non-degenerate linear transformations mixing different doublets 4> a 
but keeping invariant the kinetic term (which includes interaction of the Higgs fields with 
the gauge sector of the model). Alternatively, they can be defined as norm-preserving 
transformations of doublets that do not change the intradoublet structure. 

Keeping the kinetic term invariant (and preserving charge conservation) requires 
conservation of the norm which results in the reparametrization transforma- 

tion to be unitary (a Higgs-family transformation) or antiunitary (a generalized CP- 
transformation) : 



U : (j) a ^- U a b4>b or U C p = U ■ J : 



Pa ^ U ab (j)\ 



'b ■ 



(5.1) 



with a unitary matrix U a b- The antiunitary transformations Ucp = U ■ J, will be studied 
in Chapter [6} 

In this Chapter we focus on the unitary transformations U . A priori, such trans- 
formations form the group U(N). However, the overall phase factor multiplication is 
already taken into account by the U(l)y from the gauge group. This leaves us with the 
SU(N) group of reparametrization transformations. 



U(N) 



SU{N) 



V 



■ e 



det=l 



7 



U(l) 

i0o 



(5.2) 



where /3 = A = on - f3 and S/3 n = 0. 

Then group SU(N) has a non-trivial center Z(SU(N)) = Z N generated by the diag- 
onal matrix exp(2ni/N) ■ 1^, where Iat is the identity matrix. Therefore, the group of 
physically distinct unitary reparametrization transformations G is 



G = PSU(N) ~ SU{N)/Z N . 



(5.3) 



5.3 Finding Abelian groups in NHDM 

In this Section we introduce the strategy to find all realizable Abelian subgroups of 
the NHDM potential for any N. This strategy can be outlined as follows: we first 
describe maximal Abelian subgroups of SU(N), then we construct the maximal Abelian 
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subgroups of PSU(N). Next we explore the realizable subgroups of this maximal Abelian 
subgroup. 

Before presenting the strategy of constructing the maximal Abelian groups in SU(N), 
we start with an example which would facilitate understanding the main idea. 



5.3.1 Heuristic example 

Consider a given potential V of the N Higgs doublet model. We want to find phase 
rotations which leave this potential invariant. These phase rotations form the group 
[[/(l)]^, the group of all diagonal unitary N xN matrices acting in the space of doublets. 
This group is a subgroup of the group of all possible phase rotations of doublets U(N), 
[U(1)] N C U(N), and can be parametrized by N parameters acj G [0, 27r): 

diag[e iQl , e ia \ e iaN ] (5.4) 

The potential V is a collection of k monomial terms each of the form {4>\4>b) or (0t0;,)(0j0d). 



Upon a generic phase rotation (5.4), each monomial term gains its own phase rotation. 
For example, with a ^ b gains the phase a\, — a a , (0l0fe)(0a0 c ) with a,b,c all 

distinct gains the phase a\, + a c — 2a a , etc. 

So, each monomial gets a phase rotation which is a linear function of a's with integer 
coefficients: Y2f=i m j a j- The vector of coefficients rrij can be brought by permutation 



and overall sign change to one of the following forms: 



(1,-1,0,...) : ((f>%) 



(2,-2, 0, ... 
(2, -1, -1, ... 

T, l, -l, -l, ... 



UbWM (5-5) 



Note that in all cases Y2j m j = 0' which simply means that each bilinear is U(l)- 
symmetric. 

The phase transformation properties of a given monomial are fully described by its 
vector rrij. The phase transformation properties of the potential V, which is a collection 
of k monomials, is characterized by k vectors mi t j, m^j, ■ ■ ■ , mkj, each m it j being of one 



of the types in (5.5). 

If we want a monomial to be invariant under a given transformation defined by 
phases {ctj}, we require that Yl!j=i m i a j = ^ 1xn with some integer n. If we want the 
entire potential to be invariant under a given phase transformation, we require this for 
each individual monomial. In other words, we require that there exist k integers {rii} 
such that the phases {ctj} satisfy the following system of linear equations: 

N 

rriijaj = lurii , for all 1 < % < k (5-6) 

i=i 
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Solving this system for {a,} yields the phase rotations that leave the given potential V 
invariant. 



One class of solutions can be easily identified: if all a>j are equal, a>j = a, then (5.6) 
with rii = is satisfied for any a. These solutions form the U(l) subgroup inside 
and simply reflect the fact that the potential is constructed from bilinears (<p' a <f>b)- These 



solutions become trivial when we pass from the U (N) to the SU(N) group of transforma- 



tions. However, there can exist additional solutions of (5.6). They form a group which 
remains non-trivial once we pass from U(N) to SU(N) and further to PSU(N), which 
is the group of physically distinct Higgs-family reparametrization transformations. It is 
these solutions that we are interested in. 

In order to find these solutions, we note that a matrix with integer entries can be " di- 
agonalized" by a sequence of elementary operations on its rows or columns: permutation, 
sign change, and addition of one row (column) to another row (column). " Diagonaliza- 
tion" for a non-square matrix means that the only entries my that can remain non-zero 
are at i = j. After that, the system splits into k equations on N phases of the generic 
form 

m^a- = 2irn' i , a[ G [0, 2ir) , n { G Z (5.7) 
with non-negative integers m^. 



• If rriii = 0, this equation has a solution for any af, the 2-th equation gives a factor 
U(l) to the symmetry group of the potential. 

• If m^i = 1, then this equation has no non-trivial solution, and the i-th equation 
does not contribute to the symmetry group. 

• If m^i — di > 1, then this equation has d{ solutions which are multiples of q.{ = 
27r/di, and the i-th equation contributes the factor Z^ to the symmetry group. 



The full symmetry group of phase rotations is then constructed as direct product of these 
factors. 

Thus, the task reduces to studying which diagonal values of the matrix rriij can arise 
in a model with iV doublets. For small values of iV this task can be solved explicitly, 
while for general iV one must rely upon subtle properties of my which stem from (5.5). 



5.3.2 Maximal Abelian subgroups of PSU(N) 

In the previous Section we outlined the main idea of the strategy. However, we worked 
there in the group U(N), while the reparametrization group is PSU(N). 

We start by reminding the reader that a maximal Abelian subgroup of G (5.3) is an 
Abelian group that is not contained in a larger Abelian subgroup of G. In principal, 
there can be several maximal Abelian subgroups in G, any subgroup of G must be either 
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a maximal one, or lie inside a maximal one. Therefore, we first need to identify all 
maximal Abelian subgroups of G, and then study their realizable subgroups. 
In SU(N), all maximal Abelian subgroups are the so called maximal tori |48j . 

[U^f- 1 = 17(1) x 17(1) x • • • x U(l) (5.8) 

All such maximal tori are conjugate to each other, which means with two given 
maximal tori 7\ and there exist g e SU(N) such that g~ x T\g = T^. Therefore, 
without loss of generality one could pick up a specific maximal torus, for example, the 
one that is represented by phase rotations of individual doublets 

diag[e iQl , e ia \ e iaN ~\ e~ l ^ a *} (5.9) 

and study its subgroups. 

The analysis then proceeds as we explained in the previous Section, with an additional 
condition that all a's sum to zero, (Yl a N = 0). 

However, the group of distinct reparametrization transformations is PSU(N), which 
has a richer structure. It is proved that there are two sorts of maximal Abelian groups 
in PSU(N) 



Maximal tori, which will be constructed in Section 5.3.3 



• Certain finite Abelian groups which are not subgroups of maximal tori and should 
be treated separately. 

5.3.3 Maximal tori in PSU{N) 

Here we explicitly construct the maximal torus in PSU(N). We first introduce some 
convenient notation. A diagonal unitary matrix acting in the space of Higgs doublets 



and performing phase rotations of individual doublets, such as (5.9), will be written as 
a vector of phases: 

(ax, q 2 , ■ ■ ■ , a/v-i, - ^ aj) (5.10) 

In addition, if M\, . . . , are subsets of a group G, then (M%, . . . , Mj.) denotes the 
subgroup generated by Mi U . . . U Mj, of G. Then we construct a maximal torus in 
SU(N) 

T = (U(l) 1 ,U(lh,--- ,U(l)*-i) 

where 

C7(l)i = «i(-l, 1, 0, 0, 0) 
U(l) 2 = « 2 (-2, 1, 1, 0, 0) 
C/(l) 3 = a 3 (-3, 1, 1, 1, 0) 

U(1) N ^ = a N ^(-N+l, 1, 1, 1, 1) (5.11) 
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with all at G [0, 2ir). 

For every j we have (C/(l)i, . . . , Z7 (l)^_i) fl Z7 (l)j = {e}. Therefore: 

T = Z7(l)! x U(l) 2 x ■ ■ ■ x E/(l)iv_i (5.12) 

In particular, any element u G SU(N) can be uniquely written as 

u = ui(a 1 )u 2 (a 2 ) • • • Uiv_i(aiv-i) , G (5.13) 

Moreover, the center Z(SU(N)) is contained in the last group and is generated by 
ajv_x = 2n/N. One can therefore introduce 



17(1)^ = U{1) N - X /Z{SU(N)) (5-14) 
which can be parametrized as 

( N~ l l l \ y v 

Uil^^a^ — (5.15) 
where ojtv-i G [0, 2ir). Therefore, the maximal torus in PSU(N) is written as 



T = U(l) 1 xU(l) 2 x---xU(l) N _ 1 (5.16) 
5.3.4 Identifying the symmetries of the potential 



Next we study which subgroups of the maximal torus T in (5.16) can be realizable in 
the scalar sector of NHDM. 

We start from the most general T-symmetric potential: 

V(T) = -J^mKfya) + J2 X ^^a)(4^) +J2 X MM(4^a) (5.17) 

a a,b a^b 

Each term in this potential transforms trivially under the entire T. The important 
fact is that a sufficiently general potential of this form has no other unitary symmetry. 



In fact, when we start from the T-symmetric potential (5.17) and add more terms, we 
will never generate any new unitary symmetry that was not already present in T |48j. 
This is the crucial step in making sure that the groups described below are realizable. 
Our task now is to find which subgroups of T can be obtained in this way. 



Consider a bilinear <p\<pb where a ^ b. It gains a phase change under T (5.16) which 
linearly depends on the angles af 

4>l(j) b -> exp[z(p af ,ai + q ab a 2 H h t ab a N -i)) ■ 4>l4> b (5.18) 

with some integer coefficients p ab , q a b, ■ ■ ■ , tab- Note that all coefficients are antisym- 
metric in their indices: p ba = —p a b- These coefficients can be represented by real anti- 
symmetric matrices with integer values, or graphically, as labels of the edges of iV — 1 
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oriented graphs, one for each U(l) group. Each such graph has iV vertices, correspond- 
ing to doublets <f) a ; all vertices are joined with arbitrarily oriented edged, orientation 
indicated by an arrow. An edge oriented from (fi b towards 4> a (edge b — > a) is associated 
with the bilinear <p\(j)b and is labelled by pt, a in the first graph, q ba in the second graph, 



etc. Examples of these graphs are shown for 3HDM and 4HDM in Section 5.5 

The Higgs potential is a sum of monomial terms which are linear or quadratic in 
<$ a $b- Consider one such term and calculate its coefficients p, . . . ,t. Let us first focus on 
how this term depends on any single U(l)i subgroup of T. There are two possibilities 
depending on the value of the i-ih coefficient: 

• If the coefficient k in front of oti is zero, this terms is [/(l)j-symmetric. 

• If the coefficient k 7^ 0, then this term is symmetric under the 1i\k\ subgroup of 
[/(l)j-group generated by phase rotations by 2ir/\k\. 



However, even if a given monomial happens to have finite symmetry groups with respect 
to each single U(l)i, its symmetry group under the entire T is still continuous: for any 

set of coefficients one can adjust angles aij in such a way that p a b&\ + V tab&N-i — 0. 

Therefore, to make sure the potential is symmetric under a finite symmetry group, we 
must sum at least N — 1 such monomials. 

When studying symmetries of a given term or a sum of terms, we cannot limit 
ourselves to individual U(l)t groups but must consider the full maximal torus. 

The strategy presented below guarantees that we find all possible realizable subgroups 
of the maximal torus T, both finite and infinite. 

Consider a Higgs potential V which, in addition to the T-symmetric part (5.17) 
contains k > 1 additional terms, with coefficients pi, qi, ...t\ to pk, qk 
potential defines the following (N — 1) x k matrix of coefficients: 



...U 



This 



X(V) 



f Pi qi ■■■ h\ ( ?»u 

P2 <?2 ••• *2 



\ Pk qk 



t 



k J 



rrii 



.N-l \ 



\ m k ,i 



(5.19) 



mk,N-i J 



Here the second form of the matrix agrees with the notation of (5.6). The symmetry 
group of this potential can be derived from the set of non-trivial solutions for on of the 
following equations: 



X(V) 



' aii i 




( 2-nni \ 


\ ajv-i J 




\ 2vrn fc J 



(5.20) 



There are two major possibilities depending on the rank of this matrix. 
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Finite symmetry group 

If rankX(V) = N — 1, then there is no non-trivial solution of the equation (5.20) 
with the trivial right-hand side (i.e. all n.j = 0). Instead, there exists a unique 
solution for any non-trivial set of and all such solutions form the finite group 
of phase rotations of the given potential. 

To find which symmetry groups can be obtained in this way, we take exactly 
N — 1 monomials, so that the matrix X(V) becomes a square matrix with a non- 
zero determinant. After diagonalizing the matrix X(V) with integer entries [by 
swapping rows (columns), adding and subtracting rows (columns)], the matrix 
X(V) becomes diag(di, . . . , (ijv-i), where di are positive non-zero integers. This 



matrix still defines the equation (5.20) for a' { , which are linear combinations of 
ojj's and with n! i G Z. Therefore, the finite symmetry group of this matrix is 
Z f j 1 x • ■ • x Zd iV _ 1 (where Zi means no non-trivial symmetry). 

Note also that each of the allowed manipulations conserves the absolute value of 
the determinant of X(V). Therefore, even before diagonalization one can calculate 
the order of the finite symmetry group as | detX(V)|. 

This derivation leads us to the strategy that identifies all finite subgroups of torus 
realizable as symmetry groups of the Higgs potential in NHDM: write down all 
possible monomials with iV-doublets, consider all possible subsets with exactly N — 
1 distinct monomials, construct the matrix X for this subset and find its symmetry 
group following the above scheme. Although this strategy is far from being optimal, 
its algorithmic nature allows it to be easily implemented in a machine code. 

Continuous symmetry group 

If r&nkX(V) < N — 1, so that D = (N - 1) - vankX(V) > 0, then there exists a 



.D- dimensional subspace in the space of angles which solves the equation (5.20) 
for the trivial right-hand side. One can then focus on the orthogonal complement of 
this subspace, where no non-trivial solution of the equation is possible, and repeat 
the above strategy to find the finite symmetry group Gr, in this subspace. The 
symmetry group of the potential is then [[/(1)] D x Go- 



( di 



X(V) 



di 



\ 



(5.21) 



o 7 
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5.4 Abelian symmetries in general NHDM 



The algorithm described in Section ( |5.3[ ) can be used to find all Abelian groups realizable 
as the symmetry groups of the Higgs potential for any N. We do not yet have the full 
list of finite Abelian groups for a generic N presented in a compact form, although we 
put forth a conjecture concerning this issue, see Conjecture [4] below. However, several 
strong results can be proved about the order and possible structure of finite realizable 
subgroups of the maximal torus. 

Throughout this Section we will often use n := N — 1. Also, whenever we mention 
in this Section a finite Abelian group we actually imply a finite realizable subgroup of 
the maximal torus. 



5.4.1 Upper bound on the order of finite Abelian groups 



It can be expected from the general construction that for any given N there exists an 
upper bound on the order of finite realizable subgroups of the maximal torus in NHDM. 
In this Section we prove the following theorem: 

Theorem 1. The exact upper bound on the order of the realizable finite subgroup of 
maximal torus in NHDM is 

\G\ < (5.22) 

Before presenting the proof, le t us fi rst develop some convenient tools. First, with 

nstruc 

(Pi, Qi, 



the choice of the maximal torus (|5.16|), we construct n = N — 1 bilinears 



i = 1, ... ,n. The vectors of coefficients a 
easily written: 



U) defined in (5.18) can be 



a x = 


(2,3,4,.. 


.,71,1) 


a 2 = 


(1,3,4,.. 


.,71, 1) 


a 3 = 


(1,2,4,.. 


.,71, 1) 


a n = 


(1,2,3,.. 


.,71- 



(5.23) 



One can use these vectors to construct the n x n matrix A: 

( a, \ 



A 



det A = 1 



(5.24) 



\a n J 



From the unit determinant we can also conclude that after diagonalization the matrix A 
becomes the unit matrix. 
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Now consider a bilinear {4>\4>j) with i,j ^ 1; its vector of coefficients can be rep- 
resented as a,j-i — a,_i, [(4>l4>j) = ((p\(f>i)(4>\4>j)}. More generally, for any monomial 
{4>\4>j){4>\4>rn) with any i,j,k,m, the vector of coefficients has the form Oj_i — a^-i + 
a m _i — afe_i, where ao is understood as zero. This means that the vector of coefficients 
of any monomial can be represented as a linear combination of a's with coefficients 0, 
±1 and ±2: 



U 3 f 



dj-x — di-i 

a 3-l 

2a ? _ 1 

d y x T" di—\ 



2a 



2oj_i 



->■ 
— >• 
— >■ 



+ 1,-1 
+1,0 
+2,0 
+1 

+2,-2 



(5.25) 



Since the X-matrix is constructed from n such vectors, we can represent it as 

X ik = CijA jk (5.26) 

The square n x n matrix c^- can contain rows only of the following nine types (up to 
permutation and overall sign change): 



type 1 




M(<M0i) 


(1, o, ...,o) 




type 2 


( 




(2, 0, • • • , 0) 




type 3 


(A 




(1, i,o, ...,0) 




type 4 


(A 




(i, -i,o, ...,o) 




type 5 


(4 




(2,-l,0,...,0) 




type 6 


(4 


^•)(4^) 


(i, 1,-1, o, 


0) 


type 7 


( 




(2,-2, 0,-..,0) 




type 8 


(</>l<5 




(2,-1,-1, 0, 


,0) 


type 9 


(4 




(1, 1, -1, -1, 0, 





(5.27) 



0) 

It follows from (5.26 ) and ( 5.24[ ) that det X = det c • det A = det c. Therefore, the order 
of any finite group is given by the module of determinant of c: \G\ = | det c|. 

Let us also note two properties of the strings of type 1-9. Take any such string of 
length n, X( n ) = (xi, X2, ■ ■ ■ , x n ), which is obtained from (5.27) by an arbitrary permuta- 
tion and possibly an overall sign flip. Then any substring x/ n -i) — (xi, ■ ■ ■ , %k+i, • • • , 
obtained by removing an arbitrary element Xk is also of type 1-9. Moreover, the ele- 
ment removed can be added at any place, and still the string x',^. = (xi, . . . , x m + 
Xk, • • • , Xk-i, Xk+i, • • • , x n ) remains of type 1-9. Both properties can be proved by direct 
inspection of all the strings. 

Now we are ready to prove Theorem [T] 
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Proof. We prove by induction. Suppose that for any (n — 1) x (n — 1) square matrix 
D n _i whose rows are strings of type 1-9, its determinant = det D n _i is limited 
by \d n -i\ < 2 n_1 . Take now any n x n matrix D n constructed from the same family 
of strings and compute its determinant d n by minor expansion over the first row, with 
^i-ij dn-ii • ■ ■ > being the relevant minors. The procedure then depends on what type 
the first row is. 

1. If the first row is of type 1 or 2, then \d n \ < 2|<i^ 1 2 1 | < 2™. 

2. If the first row is of type 3 or 4, then \d n \ = \d^\ ± < + \ d n-i\ < 2 "- 

3. If the first row is of type 5, then we permute the columns so that it becomes exactly 



as in (5.27) and then add the second column to the first one. Then, the first row 
becomes (1, —1, 0, • • ■ , 0). The first minor does not change and contains rows such 
as (x 2 , x$, . . . , x n ), while the second minor contains rows (xi + x%, 23, . . . , x n ). 
Due to the properties discussed above, these strings are also of type 1-9, therefore 
the induction assumption applies to both minors. We therefore conclude that 
\d n \ = \dV 1 + d%\\<2\ 

4. If the first row is of type 6, then repeat the same procedure with only one change, 
that is we add the second column to the third. The first row becomes (1, 1, 0, . . . , 0), 
while the other rows have generic form (x\, X2, X2 + 23, 24, . . . ). The first minor 
contains rows of the form (x 2 , 2 2 +2 3 , x 4 , . . . ), which is equivalent to (x 2 , 23, 2 4 , . . . ) 
while the second minor contains (21, 2 2 + 2 3 , 24, . . . ). Both rows are of type 1-9, 
therefore \d n \ = [d^ - dj^il < 2 n . 

5. If the first row is of type 7, we follow the procedure described for type 5 and get 
\d n \ = 2|4 1 2 1 | < 2 n . 

6. If the first row is of type 8, we add the second and the third columns to the 
first, so that the first row becomes (0, —1, —1, 0), while the other rows are 
of the form (x\ + X2 + X3, x 2 , 23, 24, ... ). The second minor then contains rows 
(x\ + x 2 + 23, X3, X4, . . . ) which are equivalent to (x± + x 2 , 23, 24, . . . ), being of 
an allowed type. The third minor contains (21 +x 2 + 23, x 2 , 24, . . . ) equivalent to 

(21+23, 2 2 , 2 4 , ... ), again of an allowed type. Therefore the induction assumption 

(2) (3) 

applies to both minors, and we conclude that \d n \ — \d) l _ 1 — d y n _{\ < 2 n . 

7. If the first row is of type 9, we add the third column to the first and the fourth 
column to the second. The first row turns into (0, 0, —1, —1, . . . , 0), while the 
other rows become of the form (21 + 23, 2 2 + 24, 23, 24, . . . ). The third minor is 
built of rows (21+23, 2 2 +24, 24, . . . ), which are equivalent to (21+23, 2 2 , 24, . . . ), 
while the fourth minor is built of rows (21 + 23, 2 2 + 24, 2 3 , . . . ) equivalent to 
(21, 2 2 + 2 4 , 2 3 , . . . ). Both are of the allowed type, so the induction assumption 
applies to both minors, and \d n \ — |rf n 4 2i — c^i-il — 2 n - 
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Therefore, \d n \ < 2 n follows for any type of the first row, which completes the proof. □ 
5.4.2 Cyclic groups and their products 

In this Section we prove two propositions which show that a rather broad class of finite 
Abelian groups are indeed realizable as the symmetry groups of the Higgs potential in 
the NHDM. First we show which cyclic groups are realizable and then consider direct 
products of cyclic groups. 

Before we proceed to the main results, we need to prove that whenever we have the 



decomposition (5.26), it is indeed sufficient to diagonalize the matrix c^- instead of full 
matrix X in order to reconstruct the finite symmetry group. Let us recall that the matrix 
X defines a system of linear equations on angles ctf 

Xijctj = c ik A kj aj = 2-nri (5.28) 

This system remains the same upon 

• simultaneous sign change of the fc-th column in c and fc-th row in A 

• simultaneous exchange of two columns in c and two rows in A 

• simultaneous summation of two columns in c and subtraction of two rows in A 

• sign changes of the fc-th column in A and 

• exchange of columns i and j in A and exchange -H- aj 

• summation of two columns in A and subtraction of two a's 

• sign changes of the fc-th row in c and of the integer parameters r/. 

• exchange of two rows in c and of two r's 

• summation of two rows in c and of two parameters r 

It means that all allowed manipulations with integer matrices described earlier can be 
used for A and c. 

The next step is to diagonalize A; the matrix c becomes modified by a series of 
allowed transformations. But diagonal A is equal to the unit matrix. Therefore, we are 
left only with c in the system of equations. We then diagonalize c and construct the 
symmetry group from its diagonal values. 

This means that in order to prove a given group is realizable, we simply need to give 



an example of matrix constructed from rows of type 1-9 (5.27) , which yields the 
desired diagonal values after diagonalization. 
Now we move to two propositions. 
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Proposition 2. The cyclic group Z p is realizable for any positive integer p < 2 Ti 
Proof. We start with the following n x n matrix Cif 



( 2 -1 
2-10 
2 -1 



C2« 



\ 















2 -1 
2 



(5.29) 



We could manipulate this matrix, by multiplying the last column by 2 and adding it to 
the n — l'th column, and repeating this procedure for other columns, we arrive at 








-1 





•• 


• 





\ 










-1 


•• 


• 


















-1 •■ 


• 


















•■ 


• 


-1 




V 


2" 


2 n-l 


2 n-2 


2«-3 . . 


• 2 2 


2 


/ 



Now we could vanish the elements 2 n_1 , • • • , 2 2 , 2 in the last row, by adding other rows 
to it as many times as needed. Then we swap the rows (— l) n times to bring the matrix 
into a diagonal form, which also kills the minus sign in front of all the l's; 



C-2" 



( 2 n 

1 

1 





\ 





1 



which produces the finite group Z^n 



Now, let us modify C2« in (5.29) by replacing the zero at the left lower corner by — 1. 



— 1, then the 



Then the same transformation leads us to the group Z2«-i. 

If, instead, we replace any zero in the first column by —1, Cki 
diagonalization procedure leads us to the group Z p with p = 2 n — 2 n ~ k . 

Now, generally, take any integer < q < 2 n and write it in the binary form. This 
form uses at most n digits. Write this binary form as a vector with n components and 
subtract it from the first column of the matrix C2». Then, after diagonalization, we 
obtain the symmetry group Z p with p = 2 n — q. This completes the proof. □ 
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Just to illustrate this construction, for n 

( 1 \ 



1 



Q 



23 



10 



10111c 



1 

V 1 / 



5 and q - 
( 1 



23 we write 



c 25 - q 







V 



i o 
-i 

2 





2 









-1 
2 










2 1 -2 C 



Repeating the above procedure we obtain at the lowest left corner 2 5 — 2 4 — 2 2 
32 — 23 = 9, which gives the Z 9 group. 

Now we show that not only cyclic groups but many of their products can be realized 
as symmetry groups of some potential. 

Proposition 3. Let n = Yli=i n i be a partitioning of n into a sum of non-negative 
integers n^. Then, the finite group 



G 



is realizable for any < pi < 2 ni . 



x 



• x Z 



(5.30) 



Proof. Let us start again from the matrix (5.29). For any partitioning of n one can turn 
this matrix into a block-diagonal matrix by replacing some of the — l's by zeros. The 
matrix is then represented by the smaller square blocks of size n\, n 2 , . . . , n^. The i-th 



block has exactly the form of (5.29), with n replaced by rij. Therefore, within this block 



one can encode any cyclic group Z Pi , with < pi < 2 n \ Since each block can be treated 



independently, we can encode any group of the form (5.30) 



/ 2 -1 

2 














2 -1 


















\ 









2 




-1 

2 J 



(5.3i; 



□ 



Let us note that this proposition does not yet exhaust all possible finite Abelian 
groups with order < 2 n . For example, for n = 5, we can think of the group Z 5 x Z 5 
whose order is smaller than 2 5 = 32. However, there exists no partitioning of 5 that 
would lead to this group by applying the proposition just proved. At this moment, we 
cannot prove whether such groups are also realizable. However, we conjecture that they 
are. Together with the earlier conjecture on the upper bound on the order of the finite 
group, we can now propose the following conjecture: 
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Conjecture 4. Any finite Abelian group with order < 2 is realizable in NHDM. 

If proven true, this conjecture will give the complete classification of realizable sub- 
groups of the maximal torus in NHDM. 

5.5 Examples of 3HDM and 4HDM 

In this Section we illustrate the general strategy presented above with the particular 
cases of 3HDM and 4HDM. We give the full list of Abelian groups realizable as symmetry 
groups of the Higgs potential and show explicit examples of such potentials. 

5.5.1 Unitary Abelian symmetries of the 3HDM potential 

In 3HDM the representative maximal torus T C PSU (3) is parametrized as 

T = U(l) 1 xU(l) 2 , I7(l)i = a(-l, 1, 0) , U(l) 2 = f3 (~, (5.32) 

where a, (3 e [0, 2tt). The most general Higgs potential symmetric under T is 

V = -mli^-mli^-ml^lfo) (5.33) 
+An(0l0i) 2 + A 22 (40 2 ) 2 + A 33 (0 3 3 ) 2 

+Al 2 (0 t 10l)(0 2 >2) + A 23 (402) (4^3) + Ai 3 (0 3 3 )(0 t 10l) 

+\' 12 \44>2\ 2 + a 23 |0 2 3 | 2 + K 3 \4<Pi\ 2 

There are six bilinear combinations of doublets transforming non-trivially under T with 
coefficients 





V 


Q 


(0200 


-2 


-1 


(0302) 


1 





(0l03) 


1 


1 



(5.34) 

and their conjugates whose coefficients p and q carry opposite signs with respect to 



(5.34). These coefficients are shown graphically as labels of the edges of two oriented 



graphs shown in Figure (5.2). 



Before applying the general strategy described in the previous Section, let us check 
which finite symmetry groups can arise as subgroups of either U(l)i or £7(1) 2 individually. 
Consider first 2/ The coefficient p of any monomial can be obtained either by 



picking the labels from the first graph in Figure (5.2) directly or by summing two such 
labels, multiplied by —1 when needed. In this way we can obtain any \p\ from to 4. 
For a monomial with p = 0, the symmetry group is the entire U(l)i. For \p\ = 1, there is 
no non-trivial symmetry. For \p\ = 2,3,4, we obtain the cyclic group Zui. In each case 
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Figure 5.2: Coefficients p and q as labels of triangles representing bilinears in 3HDM 

it is straightforward to construct the monomials with a given symmetry; for example, 
((/>i</> 3 )(0203) and its conjugate are [/(l)i-symmetric, while (4i<p2)(4ih), ((^Mttlfc) 
and their conjugates are symmetric under the Z 3 -group with a generator 

a = (to, 1, w -1 ) , to = exp (2ni/3) (5.35) 

In the case of the group U(l) 2 , the labels can sum up to \q\ = 0, 1, 2. Thus, the largest 
finite group here is Z 2 . 

As we mentioned in the previous Section, one cannot limit oneself to subgroups of 
individual U(l)i factors or to direct products of such subgroups. In order to find all 
realizable groups, one has to write the full list of possible monomials and then calculate 
the symmetry group of all distinct pairs (for iV = 3) of monomials. For example, if 
the two monomials are v\ = Ai(0{02)(<M03) an d v% = A2(020i)(02 < / , 3)' then the matrix 
X(vi + v 2) has form 

*<«. + *)= i,)-^ J) (5.36) 

and it produces the symmetry group Z3. The solution of the equation 

+ = (5-37) 



P J \ 2vrn 2 

yields a = 2n/3 • k, /3 — 0, which implies the transformation matrix of the doublets 



(5.35). 



In 3HDM there are, up to complex conjugation, three bilinears and nine products of 
two bilinears transforming non-trivially under T. 

The full list of monomials in this case is as follows: 

• Linear: , {4M, (4ifo) 

. Quadratic: (4 2 <Pi) 2 , (40 2 ) 2 , (<p\fo)\ (4<f>i)(4fa), (4hMh), (4ifoM<Pi), 
(4^(4^), (4hM<Pi), (4M(4M and also (4<h)i<t>l<l>i), 
(4i4>3)(4i < fii) an d their complex conjugates 
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It is a straightforward exercise to check all possible pairs of monomials (67 pairs); in 
fact, by using permutations of the doublets the number of truly distinct cases is rather 
small. This procedure reveals just one additional finite group Z 2 x Z 2 , which arises when 
at least two terms among (0{0 2 ) 2 , (0 2 03) 2 , (030i) 2 are present. This group is simply the 
group of sign flips of individual doublets. 

Thus, we arrive at the full list of subgroups of the maximal torus realizable as the 
symmetry groups of the Higgs potential in 3HDM: 



Z 2 , Z 3 , Z 4 , Z 2 x Z 2 , £7(1), £7(1) x Z 2 , £7(1) x £7(1) 



(5.38) 



Most of these groups were identified in [49] in the search of "simple" symmetries of 
3HDM scalar potential. In that work a symmetry was characterized not by its group, as 
in our work, but by providing a single unitary transformation S and then reconstructing 
the pattern in the Higgs potential which arises after requiring that it is S-symmetric. 
In certain cases the authors of |49j found that the potential is symmetric under a larger 
group than (S), in accordance with the notion of realizable symmetry discussed above. 
The explicit correspondence between the seven symmetries Si, . . . , SV of [49J and the 



list (5.38) is the following: 



Si -> Z 2 , S 2 £7(1) realized as £7(1) 2 , S 3 ->• Z 3 (5.39) 
S 4 ->■ Z 4 , S 5 -> £7(1) realized as £7(l)i , S 6 -> £7(1) x Z 2 , S 7 ->■ £7(1) x £7(1) 



In addition to these symmetries, our list (5.38) contains one more group Z 2 x Z 2 , which 



was not found in jl9] because it does not correspond to a "simple' symmetry. 



Let us explicitly write the potentials which are symmetric under each group in (5.38). 



U(l) x U(l) = T 

The most general T-symmetric potential of 3HDM is given by (5.33). Every term 
in this potential is of the form (4>\^>i) or (0j0j)(0j0j) from which the coefficients of p 
and q both become zero, therefore this potential is symmetric under £7(l)i x £7(1) 2 . 

U(l) 

This group can be realized in two non-equivalent ways, which are conjugate either 



to £7(l)i or £7(1) 2 in (5.32). 



The general £7(l)i-invariant potential contains the following terms, in addition to 

Al32 3 (0 t 103)(0^3) + ^C (5-40) 



(5.33): 



while the general £7(l) 2 -invariant potential contains 



+^2323 



+ 



All23(0|0l) + A 2223 (0 2 2 ) + A 3323 ( 
_|_ \ .JA^A. \(A$. 



403) 



\2113 



(5.41) 
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It must be stressed that these potentials are written for the specific convention of 



groups Z7 (l)i and £7(1) 2 used in (5.32). In other bases, the explicit terms sym- 



metric under a U(l)i or £/(l) 2 -type groups will look differently. For example, the 
term (0 2 0i)(030i) is symmetric under a £7(l)i-type transformation with phases 
(0,a, 



-a) 



U(l) x Z 2 

Looking back at (5.40), one might be tempted to think that the true symmetry 



group of this term is not £7(1) but £7(1) x Z 2 , because this term is also invariant 
under sign flip of ^3. However, inside SU (3), a transformation with phases (0, 0, it) 
is equivalent to (— 7r, 7r, 0) which is already included in U(l)x. 

A potential whose true symmetry group is £7(1) x Z 2 is given by, in addition to 



(5.33), 



A 



2523(403) 



+ h.c. 



(5.42) 



This term is symmetric not only under the full £7(1)2, but also under (— 7r,7r, 0), 
which generates a Z 2 subgroups inside C7 (l)i- 



The potential symmetric under Z 4 contains, in addition to (5.33), the following 
terms: 

Ai3 23 (0 t i03)(403) + A 1212 (0 t 1 2 ) 2 + h.c. (5.43) 



As discussed before, the potential must contain at least iV — 1 = 2 distinct terms 
to be symmetric under a finite symmetry group. Again, this set of terms is specific 
for the particular choice of the (£7(l)i, £7(l) 2 )-basis on the maximal torus. This 
group is generated simply by a 4 = (i, —i, 1). 

Z 3 

The potential symmetric under Z 3 contains 

Al 232 (0 t 102)(0302) + A 23 13 (403X0103) + A 3 12l(030l) (020l) + h.C. (5.44) 

In fact, any pair of the three terms is already sufficient to define a Z 3 -symmetric 
potential. Note also that different (£7(l)i, f/(l) 2 )-bases on the maximal torus lead 
to the same Z3-group, because inside PSU(3) the following groups are isomorphic: 

((l,co,co 2 )) ~ ((u,u 2 ,l)) ~ ((w 2 ,l,w)), u = exp(27ri/3) . (5.45) 



Z 2 x Z 2 

This group can be realized simply as a group of independent sign flips of the 
three doublets (up to the overall sign flip), generated by a 2 = (—1,-1,1) and 



79 



^1,-1,-1) . Every term in the potential must contain each doublet in pairs. 



In addition to (5.33), the Z 2 x Z 2 -symmetric potential can contain 



A 



12121 



+ A 



23231 



+ A 



313H 



+ h.c. 



(5.46) 



Z 2 

This group can be realized, for example, as a group of sign flips of the third 
doublet, generated by a 2 = (—1, —1, 1). The potential can contain any term where 
(f) 3 appears in pairs. 



In Figure (5.3) we show a pictorial description of subgroups of unitary realizable 



Abelian symmetries in 3HDM. 




Figure 5.3: Tree of subgroup of unitary Abelian symmetries in 3HDM 



5.5.2 The Z 3 x Z 3 -group 

The only finite Abelian group that is not contained in any maximal torus in PSU(3) is 
Z3 x Z3. Although there are many such groups inside PSU(3), all of them are conjugate 
to each other. Thus, only one representative case can be considered, and we choose the 
following two generators 

1 n \ / n 1 n \ 

u = exp (5.47) 








^ 




^0 


1 






CO 





V -I 










V 





^ I 









/ 
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A generic potential that stays invariant under this group of transformations is 



V 



—m 
+Ai 
+A 3 



+ 



+ 



^3) 



+ A 
+ A 2 



+ 



> 2 92, 



+ 



I^ 2 | 2 + l40 3 | 2 + i40i 12 



+ h.c. 



(5.48) 



with real M , A , Ai, A 2 and complex A 3 . 

The potential ( 5.48[ ) is not symmetric under any continuous Higgs-family transforma- 
tion, which can be proved, for example, using the adjoint representation of the bilinears 
described in Chapter [2j However it is symmetric under the exchange of any two dou- 
blets, e.g. 0i -H- 02, generating group Z 2 . But this Z 2 group doesn't commute with 



Z 3 x Z 3 group. Therefore, the potential (5.48) is symmetric under a larger group which 



is non-Abelian. We can conclude that the symmetry group Z 3 x Z 3 is not realizable for 
3HDM. 



5.5.3 Unitary Abelian symmetries of the 4HDM potential 

In the case of four Higgs doublets the representative maximal torus in PSU(4) is T 
U(l)t x U(l) 2 x U(l) 3 , where 



^(l)i = 1, 0, 0) , U{1) 2 = 0{-2, 1, 1, 0) , Z7(l) 3 = 7 



3 111 

'4' 4' 4' 4 



(5.49) 



The phase rotations of a generic bilinear combination of doublets under T is characterized 
by three integers p, q, r; 



b\(j) b ) -> exp[i(p ab a + q ab (3 + r ab ^)}(4>\(, 



(5.50) 



There are twelve bilinear combinations of doublets transforming non-trivially under T 
with coefficients: 



(5.51) 





p 


q 


r 






-2 


-3 


-1 


(<pU 




1 








{<&< 







1 









-1 


-2 


-1 






1 


1 





(A 


* 3 ) 


1 


3 


1 



and their conjugates whose coefficients p, q and r carry opposite signs. 

These coefficients can again be represented graphically as labels of edges of three 
simplices shown in Figure (5.4). 

With four doublets, one can construct (up to conjugation) 42 monomials transforming 
non-trivially under T. The full list of monomials in this case is as follows: 
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03 1 02 03 02 03 02 

Figure 5.4: Coefficients p, q and r as labels of simplices representing bilinears in 4HDM 

• Linear: (0l0 2 ), (020s), (0304), (0l0i), (0204), (0 3 0i) 

. Quadratic: (0l0 2 ) 2 , (0 2 3 ) 2 , (0 3 4 ) 2 , (0l0i) 2 , (0 2 4 ) 2 , (^0i) 2 , 

(0|0 2 )(0^0 3 ), (0l0 2 )(0t04), (0l0 2 )(0i0i), (0l0 2 )(0 2 04), (0l02)(0t0i), 

(0 2 3 )(0t0 4 ), (0 2 3 )(0l0l), (0 2 3 )(0^04), (0|0 3 )(0|0l), (0t0 4 )(0l0l), 
(0^04X0204), (0304)(030l), (040l)(0203), (0l0l) (03>l) > (0201 ) (03>l) , 

along with their complex conjugates, and also (0{0 2 )(0|0j), (0 2 3 )(0;0i), 
(0 3 04)(0-0i), (0l0i)(0i0i), (0204) (0l0i), (0 3 0i) (0i0i) and their complex conjugates 

When we search for realizable finite groups, we pick up all possible combinations of 
three distinct monomials vi, t> 2 , t> 3 , construct the 3x3 matrix X{v\ + 1> 2 + t> 3 ), check 
that it is non-degenerate and then diagonalize it to obtain diag(c? 1 , <i 2 , d 3 ). The symmetry 
group is then Z dl x Z d2 x Z d3 . Although this brute force algorithm can give the full list 
of all realizable finite subgroups of the maximal torus in 4HDM, we can in fact apply 
the results of the Section IS~4l to find this list: 

Z fe with k = 2, . . . , 8, Z 2 x Z fc with k = 2, 3, 4, Z 2 x Z 2 x Z 2 (5.52) 

In short, these are all finite Abelian groups of order < 8. Perhaps, the most surprising 
of these groups is Z 7 , because it does not appear as a realizable subgroup of any single 



U(l)i factor [no two labels on Figure (5.4) sum up to 7]. Realizable continuous Abelian 



symmetry groups can also be found following the general strategy. 



Let us explicitly write the potentials which are symmetric under each group in (5.52). 

• U(l) x U(l) x U(l) = T 

The most general T-symmetric potential of 4HDM is: 

V = -m 2 (0|0!) - m^(0^0 2 ) - m^(0 3 3 ) - m^(0^0 4 ) 

+ A 11 (0 t 1 1 ) 2 + A 22 (0+0 2 ) 2 + A 33 (0^0 3 ) 2 + A 44 (0l04) 2 

+ A 12 (0 t 1 1 )(0 t 2 2 ) + A 13 (0 t 1 1 )(0 3 3 ) + A 14 (0 t 1 1 )(0l0 4 ) 

+ A 23 (0+0 2 )(0 3 3 ) + A 24 (0 t 2 2 )(0i0 4 ) + A 34 (0 3 3 )(0l0 4 ) 



+ A' 12 |0 t 1 2 | 2 + A' 13 |0 t 1 3 | 2 + A' 14 |0 t 1 c 



|2 



'4 1 

+ A 23 |0 2 3 | 2 + A 24 |0 2 4 | 2 + A 34 |0 3 4 | 2 (5.53) 
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which are the terms of the form (0j0i) or (0-0j)(0j0j). 



U(l) 

The general [/(l)i-invarinat potential contains, in addition to (5.53), the following 
terms: 

-774} (0^3) + Aii 43 (0|0i) + A 224 3(0202) + A3343 (0303) (0403) 

+A 4 343 (04 03) 2 + A 4 142(0i0l)(0l02) + h.C. + A 4 l 13 (040l) (0103) + h.C. 
+A 4 132(0i0l)(0 3 02) + /i.C. (5.54) 



The general [/(l^-invarinat potential contains, in addition to (5.53), the following 
terms: 



-^32(03 02) + 



All32(010l) + A 22 32 (02 02) + A 3 332(0 3 03) (0302) (5.55) 
+A 3 232 (0 3 02) 2 + A 43 24 (0l03)(0204) + h.C. + A 2 ll 3 (020l) (0103) + h.C. 



The general [/(l)3-invarinat potential contains, in addition to (5.53), the following 
terms: 



-™32(03 02) + 



All32(010l) + A 2 232 (02 02) + A 3 332 (0^0 3 ) (0302) 



+A 3 232 (0 3 02) 2 - ^43(0^3) + 
+A 4 343 (04 03) 2 - ^42(0402) + 



All 43 (0l0l) + A 22 43 (02 02) + A 33 43(0 3 03) (0403) 
All42(010l) + A 2 242 (02 02) + A 33 42(0303) (0402) 



+A 4 242 (0 4 02) 2 + A 4 ll 3 (040l)(0103) + h.C. + A1241 (01 02) (0 4 0l) + h.C. 
+A 2 113(020l)(0 t 103) +h.c. 



(5.56) 



One could write terms that are Z 2 -symmetric according to the labels of the edges 



in Figure (5.4). However, to make sure that a potential is only symmetric under a 



Z 2 group, we use the strategy in Section (5.4). We start with the following 3x3 
matrix: 

'2-1 
2 -1 ) (5.57) 
2 



C 2 3 



For any group Z„ we need to write q = 2 



N-l 



p in a binary form, which uses 



at most n = N — 1 digits. Now we write this binary form as a vector with n 



components and subtract it form the first column of the matrix (5.57). 



no 



110, 




c 2 s - q 




(5.56 
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The 3x3 c-matrix (5.58) represents the coefficients of the A-matrix (5.24). There- 
fore the matrix (5.58) could be written as: 

/ 



ai — a>2 
—ai + 2a 2 — a 3 
V 2a 3 

«2 = (0102), and a 3 = 



Recalling that a\ - 
to the elements of the matrix (]5.59|) follows: 



(5.59) 

^02), the terms corresponding 



ai — d 2 
-ai + 2a 2 — a 3 
2a 3 



(0I0 2 )(0 3 1 
(0 2 3 )(0l0 3 

(0i0 4 ) 2 



(5.60) 



Therefore a Z 2 -invarinat potential contains, in addition to (5.53), the following 
terms: 

Al23l(<M0 2 )(03>l) + A 2343 (02>3)(0k>3) + A 14 l 4 (0 t l04)(0 t l04) + /l.C (5.61) 

To find the phases of the transformation with which this Z 2 group is generated, we 
have to reproduce the X-matrix form the c-matrix (5.58). 

1 



1 



X, 



ik 



— >■ 



X 







2 




1 





f «1 \ 








^ «3 / 



(5.62) 



The v4-matrix for 4HDM (n = 3) is known from (5.23). Therefore: 

-1 



X 




2 






-1 





(5.63) 



Equation (5.20) in this case is reduced to the form: 





Solving this set of equations for a, (5 and 7 results in a = (3 



(5.64) 
and 7 = 7r • k. 



Adding the phases in (5.49), with the known values of a, (5 and 7, one finds the 

(5.65) 



phases of the Z 2 group generator: 

«2 



7T 

t 



-3,1,1,1 



One can check that a 2 lies in the center Z(SU(4)), and that the terms (5.61) are 
invarinat under a. 
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Using the same strategy used for the Z 2 group, one can write a Z 3 -invarinat po- 



tential with the following terms, in addition to (5.53): 

Al23l(0l0 2 )(40l) + A 13 43(0 t 10 3 )(403) + A 24 1 4 (404)(<M04) + h.C (5.66) 



This group is generated by the transformation with phases 

a-s = ^(-5,3,3, -1) 



(5.67) 



We use the same strategy to find Z 4 -invarinat potential. In addition to (5.53), this 
potential contains the following terms: 

Ai23i(0 t i0 2 )(0ki) + Ai34 3 (0 t i03)(0 4 03) + Xuu&lfa) 2 + h.c (5.68) 
The phases of the Z4 group generator are: 



a 4 = |(-2, 1,1,0) 



(5.69) 



Following the same method, one finds a Z 5 -invarinat potential contains, in addition 



to (5.53), the following terms: 



Al232(<M0 2 )(03>2) + A 23 43(403)(04>3) + XuutylfaMlfa) + h.C. (5.70) 

This group is generated by the transformation with phases 

(5.71) 



2tt 



«5 = —(-9,7,3,-1) 



An example of a Z 6 -invarinat potential would have, in addition to (5.53), the 
following terms: 



Al23 2 (0l0 2 )(0302) + A 2 343(0203)(0403) + Ai 4 l 4 (0l04) 2 + h.C. 



The phases of the Z§ group generator are: 

7T 



a 6 = -(-3,2,1,0) 



(5.72) 



(5.73) 
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z 7 

An example of terms that possess this symmetry group is 

<h)$M + Al34 3 (0 t l03)(40 3 ) + A 24 14(404)(0!04) + h.C. 



(5.74) 



including (|5.53|). This group is generated by the transformation with phases 

-21,19,3,-1) (5.75) 



7T 



14 



The matrix (5.57) is already the appropriate C-matrix for a Z 8 symmetry. From 



this matrix one could write the terms for a potential with this symmetry: 

Al232(0 t l0 2 )(0 3 02) + Ai34 3 (0 t 103)(0l03) + \u1a{4>[<P if + h.C. (5.76) 



in addition to (5.53). The phases of the Z$ group generator are: 

7T 



«8 



16 



;-15,21,-7,l) 



(5.77) 



Z 2 x Z 2 

To find a Z 2 x Z 2 -symmetric potential, we use the strategy described in preposition 
(|3|. Starting from the matrix (5.57), we break it into two block diagonal matrices 
by replacing the element ci 2 by zero. 



/ 2 10 



\ 







2 -1 
2/ 



The first matrix represents a Z 2 group. To get a Z 2 group out of the second matrix 
one has to follow the same strategy as before; writing 2(= 2 2 — 2) in a binary code, 
and subtract the vector of this binary code from the first column of the second 
matrix. Therefore having: 



-10 



10, 



c 2 2 - q 



1 -1 
2 



Therefore the full c-matrix that represents a Z 2 x Z 2 symmetry has the following 
form: 

'20 

1 -1 ] (5.78) 
2 
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From this matrix, one could easily derive the terms of a Z 2 x ^-symmetric poten- 
tial: 

Am2(0l0 2 ) 2 + A 13 4i(0l03)(0ki) + A 1414 (^4) 2 + fc.c. (5.79) 



in addition to (5.53). Equation (5.20) in this case has two sets of solutions for a, (3 



and 7. Therefore the Z 2 X Z2 group is generated by the two transformations with 
phases 

7T , , , 7T , _ 

(5.80) 



«22 = ^(- 3 > !> !> !) 



7T 



a 22 = ^(—5, —1, 7, — 1) 



Z 2 x Z 3 

Following the same strategy, we find the terms in a Z 2 X ^-symmetric potential 
are: 

Al21 2 (0 t i0 2 ) 2 + Ai 343 (0 t l03)(0l03) + A 3 414(404)(0i04) + fc-C (5.81) 



in addition to (|5.53|). To find the generators of this group we follow the same 

in two gen 

-3,2,1,0) 



method as before, resulting in two generators 

2tt 



a 2 3 



7T 



' a 23 — a I 5, 3, 3, 1 ) 



However 2 and 3 are mutually prime numbers, therefore Z 2 x Z 3 
group only has one generator. Adding up 023 and a 23 one gets; 



(5.82) 



Z(t, so this 



4 = ^(9,-7,1,-3) 



(5.83) 



which is equivalent to the generator of Z Q upon phase rotations. 



Z 2 x Z 4 

Using the same method one could write the terms of a Z2 x ^-symmetric potential: 



Al212(0l0 2 ) 2 + Ai 343 (0l03)(0l03) + Ai 4 i 4 (0 T 1 4 ) 2 + h.C. 



a, 



(5.84) 



in addition to (5.53). This group is generated by the transformation with phases 

71 , . 7T 



«24 = g(-7, 1,5, 1) , (I 



/ 

21 



(-5,-1,3,3) 



(5.85) 



Z 2 x Z 2 x Z 2 

To obtain this symmetry, one has to break the matrix (5.57) into three block 



diagonal matrices, by replacing all — l's with zero. Therefore the simple matrix: 



(5.86) 
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results in a x Z2 x ^-symmetric potential, with terms: 

Al212(<M</> 2 ) 2 + XlS^lh) 2 + Ai41 4 (0 t l04) 2 + h.C. 



(5.87) 



in addition to (5.53). 

The phases of the Z 2 X Z 2 X Z 2 group generators are: 

a 222 = |(-3, 1,1,1) , a 222 = |(-3,3,l,l) , ^22 = ^-1,1,0,0) (5.88) 
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Chapter 6 

Abelian groups containing 
anti-unitary transformations in 
NHDM 



In this Chapter we use the strategy to find all Abelian unitary symmetries in NHDM, 
introduced in Chapter [5j and extend it to groups that include anti-unitary (generalized 
CP) transformations. We use the Higgs potential of the general NHDM and the notion 
of realizable symmetries from Chapter [5] First we introduce the group of reparametriza- 
tion transformations for anti-unitary transformations. Then we describe the strategy to 
embed an Abelian group of unitary transformation into a larger Abelian group which 
includes anti-unitary transformations. Next as an example we apply this strategy to 



3HDM for all group found in the list 5.38 The results of this Chapter are published in 



6.1 The group of reparametrization transformations 

To find the symmetries of the potential, we focus on the reparametrization transforma- 
tions. A reparametrization transformation must be unitary (a Higgs-family transforma- 
tion) or antiunitary (a generalized CP-transformation): 

U : (p a h-> Uatfo or U C p = U ■ J : <p a ^ U ah <p\ (6.1) 

with a unitary matrix U ao . The transformation J = CP acts on doublets by complex 

conjugation and satisfies J 2 = 1. 

The unitary transformations U were studied in detail in Chapter [5j 

Consider now antiunitary transformations Ucp = U ■ J, U G U(N). One can define 

an action of J on the group U (N) by J ■ U ■ J, where 

(J . U ■ J)(j) = JU(J(f)) = JU((f)*) = J(U(f)*) = {U(j)*)* = U*(f) J ■ U ■ J = U* (6.2) 
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The asterisk denotes complex conjugation. This action leaves invariant both the overall 
phase subgroup U(l)y and the center of the SU(N) [J maps a unitary transformation 
onto another unitary transformation in the same space]. Therefore, we can consider only 
such U CP = U ■ J that U G PSU(N), where PSU(N) = SU(N)/Z n is the group of 
physically distinct unitary reparametrization transformation described in Chapter [5j 

Once this action is defined, we can represent all distinct reparametrization transfor- 
mations as a semidirect product 

G = PSU(N) x Z 2 (6.3) 
where J is the generator of Z 2 , and PSU(N) a normal subgroup of G. 



6.2 Finding antiunitary transformations 

As it was described in Chapter [5] all Abelian symmetries of NHDM are subgroups of 



PSU(N). We studied which subgroups of the maximal torus T (5.16) can be realizable. 
Having found the list of realizable subgroups of torus T, we can ask whether these groups 
can be extended to larger Abelian groups that would include not only unitary but also 
antiunitary transformations. Here we describe the strategy that solves this problem. 



We defined the action of the CP-transformation J on PSU(N) given by (6.2): J ■ U 



J = U* . This action can be restricted to the maximal torus T (5.16) and even further, 
to any subgroup A C T (indeed, if a G A, then a* & A). However, J does not commute 
with a generic element of T (meaning J ■ U ■ J ^ U), therefore the group (A, J) is not, 
in general, Abelian. 

In order to embed A in a larger Abelian group, we must find an antiunitary trans- 
formation J' = b ■ J which commutes with any element a £ A: 

J'a^J'y 1 = a ^ bJaJb' 1 = a <^> baT 1 ^ 1 = a <^> b = aba (6.4) 

The last expression here is a linear matrix equation for the matrix b at any given a. 
Since a is diagonal, (aba)^ = aubijajj, so whenever aaajj ^ 1, one must set to zero. 



If at least one matrix b satisfying (6.4) is found, all the other matrices can be con- 
structed with the help of the last form of this equation. Suppose b' = xb also satisfies 
this equation; then, 



b'a^b'- 1 

)a~ 1 (xb)~ 

,-i™„ ^ „-i 



(xb)a 1 (xb) 1 = a =^ xba l b l x 1 = a =>- xax 1 



=>■ xa = ax =>• a xa = a ax =>- a xa = x (6.5) 

Therefore, x can be any unitary matrix from PSU(N) commuting with a, that is, x 
centralizes the chosen Abelian group A. 
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6.2.1 The strategy 

The strategy for embedding a given Abelian group A C T into a larger Abelian group 
Acp which includes antiunitary transformations proceeds in five steps: 



Find one b which solves the matrix equation (6.4) for all a E A. 



• Find all x G PSU (N) which commute with each element of A, that is, construct 
the centralizer of A in PSU(N). 

• Find restrictions placed on x's by the requirement that product of two different 
antiunitary transformations (x%J') and (x 2 J') belongs to A. This procedure can 
result in several distinct groups Aqp- 

• Find which terms among the A-symmetric terms in the potential are also symmetric 
under some of Acp] drop terms which are not. 

• Check that the resulting collection of terms is not symmetric under a larger group 
of unitary transformations, that is, A is still a realizable symmetry group of this 
collection of terms. 

As an exercise, let us apply this strategy to the full maximal torus T. A generic 
a G T acting on a doublet fa generates a non-trivial phase rotation ^(^i, . . . , ajv_i) 



which can be reconstructed from (5.11). The equation (6.4) then becomes 



3 #i \ / e #i 



e 



* e*** = ^ ->• e^+^bij = b^ (6.6) 



Since ipi + ^ for any the only solution to this equation is fry = for all i and 
j. This means that there is no J' that would commute with every element of the torus 
T. Thus, T cannot be embedded into a larger Abelian group T C p that would include 
antiunitary transformations. 

6.3 Examples in 3HDM 



Applying the strategy just introduced, to 3HDM for all groups from the list (5.38), we 
obtain the following additional realizable Abelian groups: 

z* 2 , z 2 xz;, z 2 xz 2 xz;, z* 4 (6.7) 

where the asterisk indicates that the generator of the corresponding cyclic group is an 
anti-unitary transformation. 
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Before we explicitly describe embedding of each of the unitary Abelian groups (5.38) 
into a larger Abelian group that contains antiunitary transformations, let's just briefly 
explain why groups such as Zg, Zg, U(l) x Z2 do not appear in this list. 

When we search for anti-unitary transformations commuting with the selected sub- 
group of T, we can indeed construct such groups at the price of imposing certain re- 
strictions on the coefficients of the T-symmetric part (5.33). This makes the potential 
symmetric under a larger group of unitary transformations, which is non- Abelian. Ac- 
cording to our definition, this means that the groups such as Zg are not realizable, 
although they are subgroups of larger realizable non-Abelian symmetry groups. 



6.3.1 Embedding U(l)i x U(l) 2 



We proved in the general case that the maximal torus cannot be extended to a larger 
Abelian group. Let us repeat the argument in the specific case of 3HDM. At the first step 
of the strategy of embedding a unitary Abelian group into an Abelian group that can 
contain antiunitary transformations, we search for a matrix b G PSU(3) that satisfies 
aba = b for all a E T where T = Z7(l)i x U(l) 2 ( |5.32 ). A generic a e T is: 



\ 



(6i 



Therefore aba 



b has the form: 

r 2ia-4i/3/3£ 



aba 



e "021 

—ia—i/3/3jj 



31 



e -tf/3 6l2 

e 2ia+2i/3/3 b22 
e ia+2if3/3 b32 



e -ia-i/3/3^ 



13 



3 M*+2i/3/3 y 



'2:-! 



'21 



33 



(6.9) 



Each element (aba)^ is equal to by up to a phase rotation -0. 



i.r 



A. 



-2a - 4/3/3 -/3/3 -a - (3/3 
-/3/3 2a + 2(3/3 a + 2(3/3 
-a - (3/3 a + 2(3/3 2(3/3 



(6.10) 



In order for an element bij to be non-zero, the corresponding phase rotation must be 
zero or multiple of 2tc. It is impossible to find such b for arbitrary a and (3. Therefore, 
T = U(l)x x ?7(1)2 cannot be embedded into a larger Abelian group. 



6.3.2 Embedding U{1) 

There are two distinct types of U(l) groups inside T: [/(l)i-type and £/(l) 2 -type. Con- 
sider first U(l)i. 
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U{l)i type 

The phase rotations for the group U(l) 1 ((3 = 0, arbitrary a) allow b to have non-zero 
entries bi 2 , b 21 and 6 33 only. With a of the form: 



(6.11) 



The equation aba = b gives: 





( e' 2ia b n 


&12 


e~ ia bi 3 


aba = 


621 


e 2ia b 22 


e °23 






e ia b 32 


&33 






b\2 





&21 














633 



6 = & 21 (6.12) 



We are free to choose b\ 2 = h 2 \ = 1, 633 = 1 (the fact that det b = — 1 instead of 1 
is inessential; we can always multiply all the transformations by the overall —1 factor). 
Then, the transformation J' which commutes with any a G U(l)i is 



J'= 1 J, aJ' = J'a VaeC/(l)i (6.13) 




Next, we search for transformations a; G PSU(3) such that x J' also commutes with any 
element a G C/ (l)i; such transformations form the centralizer of U(l)i. 

From the equation axa~ l = x, analysed with the same technique as phase rotations, 
we find that x must be diagonal. 

/ in e~ 2ia x 12 e~ ia x 13 \ f x n 

axa ^=\ e 2i «x 21 a; 22 e ia :r 2 3 =x ->■ a; = :r 22 | (6.14) 
\ e iQ x 3 i e- iQ x 32 x 33 / \ x 33 

We also know that x G PSU(3), therefore x must be of the following form: 



e 



17! 



x= I e* 72 I (6.15) 

e -»(7i+72) 

We close the group under product, by checking that the production of x J' and x'J' stays 
inside U(l)i. 

X J> . x'J' = xb ■ J ■ x'b ■ J = xb ■ x'*b (6. 16) 

Also: 



bx'*b 





1 









i ( 


V 
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Therefore: 

( e ni \ / e -*72 

e i72 e'^ | (6.17) 

e -i(7i+72) I \ e i(7(+7a) 

So, the vector of phases has the form 

( 7i-72 > 72-71 , 7i + 7 2 - 7i-72 ) (6.18) 

Which should stay inside £7(1), therefore is of the form (—a, a, 0). By introducing £ and 
7, where 

7i + 7 2 = 7i + 72 = 2£ 
7i = 7 + £ , 72 = 7 - £ 

l[ = l' + Z , 7 2 = 7'-e (6-19) 

One could write 

/ \ / \ / e* \ 

x= e^"^ = e i7 (6.20) 

V e" 2t «) / V 1 / V e" 2 * / 

We see that x belongs to a [/(l)-type group Xg. Here £ e [0, 27r/3) is an arbitrary 
but fixed parameter specifying which Xg group we take, while 7 is the running angle 
parametrizing the elements of this group. At this stage, any choice of £ is acceptable. 
Thus, we have embedded U(l)i into an Abelian group generated by (Z7(l)i, J'^) ~ U(l) x 



X r J'= ( I ( | • b ■ J (6.2L 

e" 2i « 




/" 



Therefore 



e* \ / e i5 

4' = ( e« • J' = ( e* | • J (6.22) 
e- 2i§ / \ e~ 2 ^ 

Let us now see whether this group is realizable and which £ must be chosen. The only 
U(l)i symmetric terms in the potential are (5.40) 

A(0 t 1 3 )(0 t 2 3 ) + \*(4<t>i)(4<t>2) (6.23) 
Let's see how the sum of two terms transform under ; is the phase of A. 

\ / 0i \ / 

J"0 = I e * ) • 7 ■ ( 02 ) = ( e^0i I (6.24) 



e- 2i « / V 0, / \ e-w^t 
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Therefore, 



A(0 t 1 3 )(0 2 3 ) under J? |A|e^ e -^(4</» 2 )(40 1 ) 

A*(0^i)(0^ 2 ) under g |A|e-^e 6 ^(0 t i03)(0^ 3 ) 

->• ^ A -6£ = -^ A -> e = ^A/3 (6.25) 

Therefore the sum of these two terms are invariant under J'^ provided £ = ip\/3, where 
ipx is the phase of A. 

This prescription uniquely specifies the group C PSU{3). Recall that a gen- 
eral £7(1) symmetric potential also contains the T-symmetric terms (5.33). In order to 
guarantee that the T-symmetric terms are invariant under J", we must set 

m\ = ml , An = A 22 , A13 = A 23 , A' 13 = A 23 (6.26) 

However upon this reduction of free parameters we observe that the resulting potential is 
also symmetric under under another Z 2 group generated by the unitary transformation 

d= I e~ id I (6.27) 




But this transformation does not commute with Z7" (l)i- Therefore, the true symmetry 
group of such potential is non-Abelian. According to our definition, we conclude that 
£7(1) x Z* 2 is not realizable in 3HDM. 



U(l) 2 type 

Consider now t/(l) 2 : a = 0, arbitrary j3. 



-4W3 6n e -^/ 3 6i2 




aba = ( e- i/3 <%i e 2i ^% 22 e 2i ^ 3 b 23 | = b -> by = (6.28) 



No non-trivial matrix b G PSU(3) can satisfy aba = b. Therefore, £7(1)2 cannot be 
embedded into a larger Abelian group with antiunitary transformations. Note also that 
selecting various diagonal U (1) subgroups of T, for example by setting a + (3/3 = 0, will 
result just in another version of Z7 (l)i- 



6.3.3 Embedding U(l) x Z 2 

The unitary U(l)x7* 2 emerges only when the continuous group is of the U(l) 2 -tjpe. Since 
?7(1)2 cannot be embedded into a larger Abelian group with antiunitary transformations, 
£7(1) x Z2 cannot be embedded into a larger realizable Abelian group either. 

Therefore in 3HDM all realizable continuous Abelian groups are necessarily unitary 
and cannot contain antiunitary transformations. 
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6.3.4 Embedding Z 4 

The Z 4 symmetry group with generator a arises as a subgroup of U{1)\. Therefore, b and 



J' can be chosen as before, (6.13), but the new conditions for x should now be checked. 



The vector of phases in this case is: 



71-72 » 72 - l[ , 7i + 7 2 - 7i - 72 ) 



(6.29) 



which should be a subgroup of Z4. Therefore it has to agree with (— 0), which 
could be satisfied with 71 = |n or 71 = |n + |. Thus 2 should have one of the following 
forms: 



x 












or 







-2i£ 










gif fc+if +j£ 







-2i£ 



(6.30) 



Correspondingly, two possibilities for embedding the Z 4 group arise: Z 4 x Z 2 and Zg. 
Consider first the Z 4 x Z 2 embedding. 



Z 4 x Z 2 group 



X. ■ J' 



e -*f* 

e^ fc 
1 



^4 



The Z 2 subgroup is generated by J"; 



V 






e" 2i « 



• J' 




I /; • 7 

e~ 2 ^ 



J'J 






e" 2 ^ 



• J 



(6.31) 



(6.32) 



The Z 4 symmetric terms in the 3HDM potential are (5.43) 



Al323 (0103) (0203) + Ai3 23 (030l)(030 2 ) + A m2 (0l0 2 ) 2 + A* 212 (03,0l) 



(6.33) 



As in (6.25), the first two terms are invariant under J" for £ = ^Ai 3 2 3 /3, where ^1323 is 
the phase of A1323. 

The last two terms are automatically invariant under J"; 

Am 2 (0 t 1 2 ) 2 + At 212 (0 2 0i) 2 under ^ \\ l212 \e^^ {^f + |A 1212 | e -^i 2 (0 2 x ) 2 

where ■0Ai 2 i2 * s the ph ase of Ai 2 i 2 , and therefore they are symmetric under the full Z 4 x Z 2 . 
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However, similar to the previous case, restrictions (6.26) on the T-symmetric terms 



(5.33 ), result in a potential which is symmetric under a Z 2 group generated by the unitary 
transformation 

/ e iS \ 
d= e~ iS (6.34) 




with the condition 5 = ^\ 1212 /^ resulting from; 

A 1212 (0 t 1 2 ) 2 under d IXui^^e-^^) 2 
\* 1212 (4<Pi) 2 under d | A 1212 |e"^— e 4 ^(^0 2 ) 2 

-> ^A 1212 - 45 = -^Amu -> 8 = -y (6.35) 

However, d does not commute with Z4 phase rotations, therefore, the resulting sym- 
metry group of the potential is non-Abelian. Hence, the Z4 x Z 2 symmetry group is not 
realizable. 



Z* group 

Consider now the Zg group generated by 



J'J 



J"-b 








— «7r/4 






-2i£ 



• J 



(6.36) 



with the property (J' 



///\2 



a, the generator of the Z4. 

-^/2 



W\2 








giTT/2 Q 
1 



(6.37) 



Note that upon action of J'" one of the ^-symmetric terms (5.43) (</40 2 ) 2 changes its 
sign; 

{<j>\^f under J"' e CT (<^ 2 ) 2 = -(<M</> 2 ) 2 (6.38) 



Therefore, such a term cannot appear in the potential. This means we are left with 
only one type of Z4-symmetric term in the potential. According to the discussion in 



Section 5.3.4, rankX(K) = 1 and, therefore, the potential becomes symmetric under a 



continuous symmetry group. Thus, Zg is not realizable as a symmetry group of a 3HDM 
potential. 
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6.3.5 Embedding Z3 

The Z 3 subgroup of the maximal torus T is generated by the transformation a 

2 W3 Q Q 

e 2i7r / 3 I (6.39) 
1 




This group arises from the L r (l)i-type group, therefore the representation for J' in (6.13) 
is still valid. Therefore the vector of phases has the form: 

( 71-72 > 72 - l'i , 7i + 7 2 - 7i - 72 ) (6.40) 

And has to agree with (— ^fc, ^A;, 0). With the same method as before, we see that J" 
has the form 

gi£-wr/3 




j// = e *f«r/3 I ■ J (6.41) 

e- 2i « 



Since (J") 2 = a in (6.39), this transformation generates the group Zg. The Z 3 -symmetric 



terms in the potential are (5.44) 



A 1 (0 t 1 2 )(0 t 1 3 ) + A 2 (0 t 2 3 )(0 t 2 0i) + A 3 (0 3 0!)(0 3 2 ) + h.c. (6.42) 

with complex Aj. Let's see what conditions have to be satisfied to make these terms 
symmetric under Zg. 

/ e ^-^/ 3 \ / 0i \ / e^/ 3 2 \ 
J"0 = e *+W3 • J • 2 = e i «+^/ 3 t 1 (6.43) 

V e- 2 ^ / V ^3 / V e- 2 *4 J 

Therefore 

A 1 (0 t 1 2 )(0 t 1 3 ) under J" \Xi\e^e^-^{<f>[<h)(4<h) 
A 2 (0 t 2 3 )(050 1 ) under J" \X 2 \e^ e - 1 ^ (c^M^M 
As(0 3 0i)(0 3 2 ) under g lAg^e 6 *^)^) (6-44) 
For these terms to be Z 3 invariant, the following conditions have to be satisfied: 

|Ai| = |A 2 | and ip\ 1 + ipx 2 — 3£ + tt (6.45) 
And since ip\ 3 — ~ 3£, one can deduce that 

^Ai + 4'\2 + i'Xa = n (6-46) 
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where ipi are phases of A;. 



It is, therefore, possible to have a Zg-symmetric potential by requiring (6.45) and 



applying the usual restrictions (6.26) on the T-symmetric part of the potential (5.33) 



However, just as in the previous cases, the resulting potential becomes symmetric 



under a larger unitary symmetry group, of the form (6.34), with the conditions: 

|Ai| = | As | and ip^ — ip\ 2 = 35 (6.47) 
which come from requiring the sum of the following terms to be invariant under d: 

Ai^X^s) under d lAile^- 3 "^)^) 
A 2 (40 3 )(02>i) under d |A 2 |e^e 3l5 (0 t 1 3 )(0 t i0 2 ) 
A 3 (0 3 0i)(0302) under d |A 3 |e^ (^0 2 )(^0 1 ) (6.48) 
Again this d does not commute with Z 3 , and therefore the group Zg is not realizable. 

6.3.6 Embedding Z 2 x Z 2 

The Z 2 x Z 2 C T symmetry can be realized as a group of simultaneous sign flips of any 
pair of doublets (or alternatively of independent sign flips of any of the three doublets). If 
the simultaneous sign flip of doublets % and j are denoted as Rij, then the group contains 
{1, R12, R13, R23}, with R12R13 = R23 and any (Rij) 2 = 1. Therefore the generator of 
Z 2 x Z 2 is: 

\ 

a 2 , where er^ = ±1 (6.49) 

^3 / 

From the equation axa~ l = x, we have: 

a\x n O x a 2 X l2 (TiCr 3 Xi 3 \ / X U X 12 X 13 \ 

o 2 (TiX 21 o\x 22 a 2 a 3 x 2 3 = x 21 x 22 x 23 (6.50) 
0-30-1X31 o 3 o 2 x 32 a 3 x 33 J \ x 31 x 32 x 33 J 

Since the products Oia 2 , <Ji<j 3 and a 2 a 3 can not all be +1 at the same time, therefore 
Xij = 0, i 7^ j. So, a; is a diagonal matrix, and since it has to be unitary, it has the form: 

/ e^ 1 \ 
x = e^ 2 (6.51) 
\ e"^ 1 -^ 2 / 

This group can be embedded into Z 2 x Z 2 x Z 2 , where 7L* 2 is generated by 

/ e^ 1 

Ji* = ^ 2 I ' J ( 6 - 52 ) 
\ e-^-^ 
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Up to here, any £1, 6 can be used to construct the Z2 x Z2 x Zg group. Note that this 
group can also be written as Z 2 x 7L* 2 x ^2 or ^ x ^2 x ^2 by redefining the generators. 



The Z 2 x Z 2 -symmetric potential contains the following terms (5.46); 
Ai 2 (0 t i0 2 ) 2 + A 2 



\23l 



40 3 ) 2 + A 3 i( 



40i) 2 + /i.e. 



(6.53) 



with complex Ay. Let's see with which conditions these terms are invarinat under Z 2 x 
Z 2 x Zg-symmetric potential: 



J'. 



Mi 





-M2 








J 




(6.54) 



Therefore, 



Ai2(0i02) 2 under J 7 |Ai 2 |e 



A 



31 1 



under J 7 
under J^ 7 



|A 23 [e^-e- 2 ^- 4 ^(40 3 ) 2 



I A31 |e 



#31^1+2*6(0^)2 



(6.55) 



So, we require the following conditions: 
V> 12 - 26 + 26 = -^12 

^23 - 26 - 4£ 2 = -^23 
V>31 + 46 + 26 = -V>31 



^12 
V>23 
031 



6-6 
6 + 26 
-26-6 



(6.56) 



■031 = 0, or alternatively, when the product A12A23A31 is 



Which adds up to 0i2 + 023 
real. 

Note that in contrast to the previous cases, the full T-symmetric potential (5.33) is 
invariant under Z2 x Z2 x Z?i and no larger group. Therefore this symmetry group is 
realizable. 



6.3.7 Embedding Z 2 

The unitary Abelian symmetry group Z 2 can be implemented in a variety of ways, all of 
which are equivalent. One can take, for example, the Z2 group generated by sign flips of 
the first and second doublets, R\2- The usual CP-transformation J commutes with R12. 
An analysis similar to the previous case shows that the Z2 x Z2 group generated by R12 
and J is realizable. 

It is also possible to embed the Z 2 group into Z4. Consider the following transfor- 
mation 

/ \ 

J^= _ e « J, with (J'f = R 12 (6.57) 
V e- 2i « / 
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where £ is an arbitrary but fixed parameter. The ^-symmetric potential consists of 
the T-symmetric terms (5.33) with the usual restrictions (6.26) and the specifically Z%- 
symmetric terms: 



A 5 (0l02) 2 + A 6 (0l0 2 ) [(^M - (0 2 2 )] + A 7 (0 t 1 3 )(0 t 2 3 ) 

+A 8 (0 t 1 03) 2 + A 9 (0^3) 2 + ^.C. 

Let's study how these terms transform under J 7 ; 



e* 
-e* 
e 






-2i£ 



J 




— e" 



-2i£ J.t 



Therefore, 



A, 



A 6 (0 t 1 2 ) (4^) - (0 2 2j 



) 2 under J 7 



under J 7 



A7(0l03)(0203) under J' 
A 8 (^0 3 ) 2 under 
Ag(0 2 03) 2 under 



A 5 | e ^(-0t0 2 )2 

A 6 | e ^(-0i0 2 ) [(40 2 ) 

X 7 \e^e-^(4<l>i)(4<p2 

A 8 |e^e- 64 «(0^ 2 ) 2 

A 8 |e*e- 6 ^(-40 1 ) 2 



(6.58) 



(6.59) 



(6.60) 



>i<n, 



To keep these terms invarinat, we require: 



|A 8 | — |Ag| 
^3 + ^9 = 6£ 

2^ 7 = 6£ 



(6.61) 



Or simply ip 8 + ip 9 = 2ip 7 . By adjusting A 5 and A 6 , one can guarantee that no other 



unitary symmetry appears as a result of (6.26). Therefore, this Z* 4 group is realizable 
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Chapter 7 

Zn scalar dark matter from NHDM 



The discrepancy between the masses of large astronomical objects, determined from their 
gravitational effects, and the masses calculated from the luminous matter they contain, 
led to the hypothesis of dark matter. 

In 1932, by measuring the orbital velocities of stars in the Milky Way, Jan Oort 
was the first to find evidence for dark matter, when he calculated that the mass of the 
galactic plane must be more than the mass of the material that can be seen [50J. In 
1933, while examining the Coma galaxy cluster, Fritz Zwicky used the virial theorem to 
infer the existence of unseen matter, which he referred to as "dark matter" j5H E2]- He 
calculated the gravitational mass of the galaxies within the cluster and obtained a value 
at least 400 times greater than expected from their luminosity, which meant that most 
of the matter must be dark. 

Dark matter is a currently unknown type of matter responsible for 23 percent of 
the mass-energy content of the universe. It does not emit or absorb light (or any other 
electromagnetic radiation), and therefore cannot be seen. Dark matter particle must 
be (almost) stable on cosmological time-scales [53]. Clearly its decay lifetime has to be 
larger than the age of the Universe (~ 10 18 seconds). 

Despite the astronomical evidence for the existence of dark matter (DM) [51], there 
is still no direct experimental clue of which particle can play the role of a dark matter 
candidate. None of the Standard Model particles fully explain the anomalies observed. 
However, in the models beyond SM one can find dark matter candidates, which interest- 
ingly also address other particle physics issues such as electroweak symmetry breaking 
and small neutrino masses. 

Initially DM models contained completely inert dark matter, which would not be 
destroyed in any reaction, and would only annihilate through dd* — > X$m- Many models 
were introduced in which the stability of DM is provided by a conserved Z 2 symmetry. 

The motivation for stabilization of DM using Z 2 symmetry comes from the Standard 
Model, where stability is resulted from the gauge symmetries and particle content, rather 
than from an ad-hoc symmetry; The photon is stabilized by the exact U(l) gauge sym- 
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metry, and the electron is stable as the lightest charged particle under this gauge group. 
The lightest neutrino is stabilized by the Lorentz invariance. And the proton is stabilized 
by conservation of baryon number, which results from the SM gauge symmetries and the 
gauge charges assigned to the SM particles. 

In supersymmetric models the role of the Z2 symmetry is played by the .R-parity, while 
in more phenomenologically oriented models such as the Inert doublet model [55], or the 
minimal singlet model (56], the Z2 symmetry is imposed by hand when constructing the 
Lagrangian. In these models, all the SM particles including the SM-like Higgs boson have 
positive parity, while the dark sector particles are of negative parity. The lightest among 
these negative parity particles, which we will generically denote as d here, is stable and 
represents the dark matter candidate. 

Although Z2 symmetric models avoid the decay of d to Standard Model particles, they 
allow direct two-particle annihilation dd — >■ Xsm, where Xsm is any set of SM particles. 
This certainly changes the kinetics of dark matter evolution in the early Universe and 
its relic abundance after the freeze-out. 

Trying to avoid the direct two-particle annihilation dd — > Xsm > it is natural to 
explore groups larger than Z2 and try to stabilize dark matter by conserving their discrete 
quantum number. This idea was explored in a number of papers, in which both Abelian 
[571 EBJ El] and non- Abelian [60] finite groups were used. 

One particular class of groups used to stabilize dark matter are cyclic groups 7L V 
[5T] [5*5] . With this choice, all fields are characterized by a conserved quantum number 
which we will call the Z p charge q and which is additive modulo p. Usually one assigns 
zero Z p charge (q = 0), to all the SM fields including the SM-like Higgs boson, while 
the dark matter candidates have non-zero Z p charges. This strategy prohibits the direct 
two-particle annihilation dd —> Xsm, but opens up the possibility of a new two-particle 
process such as dd — > d*XsM, which was called semi-annihilation in [61], or even multi- 
particle versions of it, "multiple semi-annihilation" . Such processes, have different impact 
on the kinetics of the dark mater abundances in the early Universe ]62j . 

In this Chapter we demonstrate that Z p -stabilized scalar dark matter can easily arise 
in multi-Higgs-doublet models. The advantage of these models is that even with few 
doublets one can get Z p with a rather large p. This fact is useful in avoiding semi- 
annihilation process dd — > d*XsM, which is allowed in a Z 3 -symmetric model and can 
be avoided by using a Z p group with p > 3. 



The structure of the Chapter is as follows. In Section 7.1 we study Z p -symmetric 
models, particularly a 3HDM example in which dark matter candidates are stabilized 
by the Z3 symmetry group. In Section 7.2 we explore Z p -symmetric 4HDMs, and the 



possibility to avoid semi-annihilation processes with large Z p , in particular we study a 
Z 7 -symmetric potential in detail. The results of this Chapter are published in 
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7.1 Z 3 -symmetric 3HDM 



Before presenting the example of a Z 3 -stabilized dark matter in 3HDM, let's list the 
conditions that an electroweak symmetry breaking (EWSB) model with Z p -stabilized 
scalar dark matter must satisfy. 

• The entire Lagrangian and not only the Higgs potential must be Z p -symmetric. 
The simplest way to achieve this is to set the Z p charges of all the SM particles 
to zero and to require that only one Higgs doublet (the SM-like doublet) couples 
to fermions. The Z p charge of this doublet must be zero, and it does not matter 
which doublet is chosen to be SM-like due to the freedom to simultaneously shift 
the Z p charges of all the doublets. 

• The Z p symmetry must remain after EWSB. This is possible when only the SM-like 
doublet acquires a non-zero vacuum expectation value. 

• If we insist on Z p -stabilization, that is, we require that not only decays but also 
2-, 3-, . . . , (p — l)-particle annihilation to SM fields are forbidden by quantum 
numbers, then the dark matter candidates must have Z p charge q which is coprime 
with p. 

To show that multi-Higgs-doublet models can easily satisfy these conditions, we start 
with the simplest model of this kind, Z 3 -symmetric 3HDM. 

A natural way to implement the Z 3 symmetry in 3HDM is to construct a potential 
symmetric under (pi — Y 02 — ► 03 — ^ 01- However, we proved in Chapter [5] that upon an 
appropriate Higgs basis change, this transformation will turn into pure phase rotations 
of certain doublets. Therefore, we use the phase rotation representations of the cyclic 
symmetry groups 0« — > e JOi 0j. Also, to keep the notation short, we will describe any 
such transformation by providing the A-tuple of phases o;,. 

A scalar potential invariant under a certain group G of phase rotations can be written 
as a sum V = Vq + Vg, where Vo is invariant under any phase rotation, while Vq is a 
collection of extra terms which realize the chosen symmetry group. The generic phase 
rotation invariant part has the form 

i ij 

while Vg obviously depends on the group. In particular, for the group Z 3 in 3HDM we 
have 

V Zs = Al(03^l)(020l) + \2(4<f>2)(4h) + \3(4foMfo) + h.C. (7.2) 

where at least two of the coefficients Ai, A2, A 3 are non-zero (otherwise, the potential 
would have a continuous symmetry). This potential is symmetric under the phase rota- 
tions 

0i^0i , 2 ^e 2 ^V 2 , 3 ^e 4 ^V 3 (7.3) 
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The vector of phases in this case has the form 



2tt 

a = — (0, 1, 2), a 3 = l (7.4) 

This notation suggests Z 3 quantum numbers 0, 1 and 2 for the first, second and third 
doublet respectively. In fact, the assignment of these charges to the three doublets is 
completely arbitrary, and the group generated by the generator a with permuted charges 



has the same action on the potential. In (7.4) we simply chose to be the SM-like 
doublet. 

Whether this symmetry is conserved or spontaneously broken depends on the pattern 
of the vacuum expectation values. If we insist on conservation of the Z 3 symmetry, we 
require that 

(0?)=^i = ^ , <03> = <$) = O (7.5) 

The potential Vo + V% 3 can have a Z3-symmetric global minimum upon appropriate 
choice of coefficients. Below we show that we do not require any constraints on A's and 
therefore no fine-tuning is required for this point to be a minimum of the potential. 

First, we note that if (v\, 0, 0) is an extremum of Vq, then it is also an extremum of 
Vo + Vz 3 , since the extra terms from V% 3 contain 0i only linearly and quadratically and 
they do not create any new linear terms in the potential after EWSB. 

Therefore, when constructing a model, one can first build Vo with a global minimum 
at (vi, 0, 0) and then add a sufficiently weak (with no strong couplings) Vi 3 so that this 
point remains a minimum. 

Now, turning to minimization of Vq, suppose that we search for the neutral minimum 
with a generic complex v.e.v. pattern (v\, t>2, 1*3). The potential would have the form; 

Vo = -mi 2 Vi 2 + XuVi 4 + (Ay + Ay) Vi 2 Vj 2 (7.6) 

Introducing pi = \vi\ 2 > 0, we can rewrite Vq as; 

V = -m 2 pi + Xup 2 + (Xij + Ay) pipj 

= -MiPi + ^AijPiPj (7.7) 

where Mj = (m 2 , m 2 ,, 777.3), and Ay is constructed from Ay and Ay. One could introduce 
Pi = (Ajj) _1 M,-, and rewrite the potential as 

Vo = ^Ay(pi - Pi){pj - Pj) + const . (7.8) 

Positivity conditions on the potential guarantee that Ay is a positive definite matrix, 
therefore its inverse exists. 
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The allowed values of the pi, (pi > 0), populate the first octant in the three- 
dimensional euclidean space. Due to Q7.8 ), the search for the global minimum (trying 



to minimize (pi — pi)) can be reformulated as the search for the point in the first octant 
which lies closest to pi in the euclidean metric defined by Ay. Clearly, if p^ itself lies 
inside the first octant (pi, P2, p 3 > 0), then the global minimum is at pi = pi. If pi lies 
outside the first octant, then the closest point lies either on the edge of the first octant, 
or at the origin. 

The global minimum is unique by the convexity arguments; since the octant pi > 
is convex, there is a unique pi that minimizes (pi — pi). 

Having established that the required vacuum pattern is generically possible, we now 
switch to a simple version of the model. This is done only to simplify the presentation 
of the argument; if needed, the calculations can be repeated for a generic potential. 
Namely, we now choose m\ > 0, m 2 ,, m\ < and take Ay = 2X Q 5ij, so the the potential 
becomes; 



V = Vo + Vz 3 (7.9) 
= -ml^M + K|(40 2 ) + \m 2 3 \(4<t> 3 ) + A [(0l0i) 2 + (4fo) 2 + (030s) 2 
+\ 1 {4<t> 1 )(44> 1 ) + A 2 (0 t 1 2 )(0^ 2 ) + A 3 (0^ 3 )(0 t i0 3 ) + h.c. 
By construction, its global minimum is at (0°) = (v/\/2, 0, 0), where v 2 = mf/Ao. 



7.1.1 Mass matrix 

In order to find the mass matrices, we write the doublets as 

*=U(, + T +iG »))> *=(*)■ * = (*) (7 - io » 

Here h is the SM-like Higgs boson, G° and G + are the Goldstone bosons, while w^, 
and z 2 , z 3 are charged and neutral scalar bosons, respectively. 

Note that in contrast to the usual practice, we describe the neutral scalar bosons in 
the second and third doublet by complex fields rather than pairs of neutral fields. The 
reason is that, by construction, the fields corresponding to the real and imaginary pairs 
of 2's have identical masses and coupling constants. In any process that can arise in 
this model, these two fields are emitted and absorbed simultaneously, so they can be 
described by a single complex field Zi. 

Now we assign the Z3 quantum numbers according to the generator a — 27r/3(0, 1, 2). 
The SM-like Higgs boson gets Z3-charge q = 0. Fields and z 2 , from the second 
doublet, have Z^-charge q — 1, while for w% and z 3 , from the third doublet, q — 2. The 
conjugates of fields w 2 and z\ have q — 2 ( = — 1 mod 3), while for , z\ q = 1 ( = — 2 
mod 3). Note that the neutral scalars with similar quantum numbers can mix. 
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Let's rewrite the potential in terms of the new parameters; 



Vo 



+ 



[-m\v + \ v 3 ] h 



(7.H) 



m; 



3A 



0„.2 



+ [-ml + \ v 2 ] (G-G + ) + 

+ l™2l( w 2^2~) + l m 3l( W 3 W 3 

+ | {z* 2 z 2 ) + \mj \(z* 3 z 3 ) + 



m| \ v 2 
"22 



A^ 2 



(23*2) 



With the value t> 2 
,2 _ r> 2 



1 , one finds the mass spectrum; The SM-like Higgs boson has mass 



"12 » "i.„± 



m 



3b 



mi = 2mf, while the masses of the charged scalar bosons are m ± 

w 2 „ a 

and the Goldstone bosons G ± , G° obviously have mass zero. 

Neutrals with equal q can mix, which indeed happens at Ai 7^ 0. We describe the 
resulting mass eigenstates by complex fields d and D (m^ <mc), both having q — 1: 



rn 



. i^l + l^sl^ 1 /,,_ 2 i i_2U2 1 A il 2 ™4 



± 2V^' m2 ' ~~ l m iD 2 + ^2- m l 



(7.12) 



where d and D are rotations of ^ 2 and Z3; 



d = cos a Z2 + sin ae 1/3 z 3 , D = — sin ae ifl z 2 + cos a 



tan 2a 



|Ai 


2 

mf 


Ao 


m 2 . 







argAi 



(7.13) 



Note that within this model we have 



m d < ™ w f » "V < m ^ > m ™± + m » 3 ± = m d + m z> 



(7.14) 



7.1.2 Higher order terms in the Z3 symmetric potential 



The triple and quartic interaction terms arising from Vq and V% 3 specify the dynamics 
of the dark matter candidates. After using a unitary gauge transformation to eliminate 
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the Goldstone bosons, higher order terms in the Z 3 -symmetric potential are: 

V 3 = A [2vhG-G + + vh 3 + vhG°G°] (7.15) 



—^Wn G + zt H —w 9 G + z% + hztzt 

y/2 V2 



+ /i.e. 



+ — ^ [-22^3 11^ + -22-22-23] + h.C. 

v2 

V A = 2\ [(G-G + ) 2 + h 4 +^(G°G ) 2 + h 2 G-G + + G°G°G-G + h 2 G° 2 

+ {w 2 W^f + (z* 2 Z 2 f + 2w 2 W%Z 2 *Z- 1 + (W3W3) 2 + (4 Z 3) 2 + 2w 3 W t 4 Z %\ 

+\ 1 [w 3 G + G + w 2 + -^=w 3 G + hz* + -^=w 3 G + G°z 2 + -^=w 2 G + hz* 



+\h 2 z*z* 3 + l -G°hz* 2 z; - ^= W2 G°G + z; - l -G°hz* 2 z; + l -G°G-z* 2 z^ 



+\2[w 3 G w 2 w 2 +w 2 G z 2 z 3 H — i=hw 3 w 2 z 2 H — ^^2-22-23 

1 i 

+A 3 [^ G_ « ; 3 fw 3 f + "yl^s" ^3" ^3 - ~^w 2 w^G°z 3 + ^^wJC - 

1 i 

+ — hz 2 z 3 z 3 —G°z 3 z 3 z 2 ] + h.c. 

v 2 v 2 

The higher order terms in the potential in terms of d and D, regardless of the coefficients 
are of the following form; 

V 3 = X [hhh] (7.16) 
+Ai [hd*d, hd*D, hdD* , hD*D] 

+A 2 [w£w 3 d,wtw 3 D,ddd,ddD,dDD,DDD] + h.c. 
+A 3 [w 2 w£d*, w 2 w£D*, d*d*d*, d*d*D*, d*D*D*, D*D*D*] + h.c. 
V A = \ [hhhh,(w 23 w+ 3 ) 2 ,(w 2 w+ + w 3 w+)(dd*,dD*,d*D,D*D), 
ddd*d*, ddD*D*,d*d*DD, DDD*D*] 
+Ai [hhd*d, hhd*D, hhdD*, hhD*D] 

+A 2 [hw 2 w 3 d, hw 2 w 3 D, hddd, hddD, hdDD, hDDD\ + h.c. 

+A 3 [hw 2 wtd*,hw 2 w+D*,hd*d*d*,hd*d*D*,hd*D*D*,hD*D*D*] + h.c. 

The lightest particle from the second and third doublets is d, and it is stabilized against 
decaying into the SM particles by the Z 3 symmetry. From the kinematically allowed 
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terms in the third order potential, we have decays such as D — > dh. Charged Higgs 
bosons wf 3 will decay into d or d* plus SM particles. 

With the Z 3 conserved quantum number, one could also write third order terms such 
as D — y dZ, and D — y w^W' , D —y w^W + , if allowed kinematically. 

If the mass splitting between the d and D is small, then these processes involve virtual 
h, Z, etc. which then decay into the SM particles. In this aspect, D decays are similar 
to weak decays. 

In the case of symmetric dark matter, the main process leading to depletion of dark 
matter after electroweak symmetry breaking is the direct annihilation dd* —y Xsm, an d 
the semi-annihilation reaction is only a correction to this process. Still, it might be 
possible that this correction leads to a sizeable departure of the kinetics of the dark 
matter in the early Universe and affects the relic abundance at the freeze-out [52] . 

The situation becomes more interesting in models with asymmetric dark matter |64j, 
in which an asymmetry between d and d* is generated at a higher energy scale. It is 
possible for example that upon electroweak phase transition almost all d's annihilate with 
d* into the SM sector, leaving behind a certain concentration of dark matter candidates 
d. 

d can also scatter elastically dd —y dd, which can initiate semi-annihilation processes 
such as dd —y d*X$M with a subsequent annihilation of d* with a d. This possibility 
originates from the following quartic terms in the scalar potential 

—=hddd cos a sina;e~ l ^(A2 cos a + A3 sinae~ J/3 ) + h.c. (7-17) 
V2 

Depending on A's, this process can be as efficient as the direct annihilation in the usual 
annihilating dark matter models, or it can be suppressed by the small coupling constants. 

Finally, the same interaction terms also generate the triple annihilation processes 
ddd — y h — y Xsm, whose rate is, however, suppressed at small densities compared to the 
semi-annihilation. 

7.2 Avoiding semi-annihilation in 4HDM 

The presence of the hddd terms in the interaction Lagrangian in the 3HDM example, 
which were responsible for the two-particle semi-annihilation process, was due to the Z3 
symmetry group. 

It is possible to avoid such two-particle semi-annihilation processes by employing a 
Z p group with larger p in a model with four Higgs doublets. 

As proved in Chapter [5j one can encode in the 4HDM scalar sector any group Z p 
with p < 8. Note that in order to avoid continuous symmetry, one must accompany the 
phase-symmetric part of the potential Vq with at least three distinct terms transforming 
non-trivially under phase rotations. 
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The full list of these realizable symmetries, their corresponding terms and their gen- 
erators are given in Section |5.5.3 In Table we remind the reader of these symmetry 
groups together with examples of the three interaction terms and the phase rotations 
that generate the corresponding group. Note that all these groups are realizable, so that 
if the three terms in each line are written down with non-zero coefficients, there is no 
phase transformation other than multiple of the generator and the common overall phase 
shift that leaves them invariant. 



Group Interaction terms Generators 
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h), (0i0 3 ), (</ 
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f (0, 0, 0, 1) 


z 3 




i>l<h)(<A<h), (< 






f (0, 1, 1, 2) 
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f (0, 1, 1, 2) 


z 5 


(0i0 3 )(0 2 3 ) 


, (0^ 2 )(0 t 1 2 ) 
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40l)(40l) 


f (0, 1, 2, 3) 


z 6 




h), (4fo){<t>\< 


K 




T (0, 1, 2, 3) 


z 7 


(40l)(40l) 




, ( 




t(0, 2, 3, 4) 


z 8 






K 


, (<M<M 2 


f (0, 2, 3, 4) 



Table 7.1: Cyclic groups realizable as symmetry groups in the scalar sector of 4HDM 



Since the potential is invariant under the common phase shift of all doublets, one can 
freely add additional equal phases to the generators shown in the third column in Table 
7.2| and, in addition, one can permute the doublets. For example, the last line of this table 

can be replaced by {(j>\(h){<A < h) ^ (030i)(040i)> (0403) 2 , (04 -> 03 -> 02 -»■ 0i ->• 04), 
which is symmetric under the Z 8 group generated by phase rotations ^ (4, 0, 2, 3) = 
f (2, -2, 0, 1) = f (2, 6, 0, 1) = f (0, 1, 2, 6). 

The patterns of phase shifts given in this table allow for construction of various dark 
sectors with different possibilities for dark matter dynamics. Here, we show that there 
are examples in 4HDM in which semi-annihilation of dark matter candidates is also 
forbidden. For this purpose we pick a Z 7 -symmetric example. 
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7.2.1 Z7 symmetric example 

Consider the Zy-symmetric 4HDM with the potential; 

v = v + v Z7 = — ^fc^t^!) -k |^li(<^+^> 2 ) -k |^ii(<^+^ 3 ) -h |^2l(<^>i<^ 4 ) 

+A 



(7.18) 



(0l0i) 2 + (40 2 ) 2 + (40s) 2 + (0l0 4 ) 2 ] 
+A 1 (0l0 1 )(^0i) + A 2 (0i0 2 )(0l0 2 ) + A 3 (0i0 3 )(403) + h.c. 
As before, we assume that only the first doublet couples to fermions (v/y/2, 0, 0, 0), 



then we expand the doublets; 



,02 



A A 



(7.19) 



^(v + h + iG°) j - y z 2 J \z 3 j" y z 4 

Similar to the previous example, h is the SM-like Higgs boson, G° and G + are the 
Goldstone bosons, while w 2 , w 3 , w 4 and z 2 ,z 3 , z 4 are charged and neutral scalar bosons, 
respectively. 

Now we assign the Z7 quantum numbers according to the generator a = 27r/7(0, 2, 3, 4). 
The SM-like Higgs boson gets Zy-charge q = 0. Fields w 2 and z 2 , from the second dou- 
blet, have Z 7 -charge q = 2, while for w 3 and z 3 , from the third doublet, q = 3, and for 
w 4 and 24, from the forth doublet, q = 4. The conjugates of fields-u^ and z\ have q = 5 
( = —2 mod 7), and and z\ have g = 4 ( = — 3 mod 7), while for w 4 , z\ q = 3 ( 
= —4 mod 7). Note that the neutral scalars with similar quantum numbers can mix. 

Let's rewrite the potential in terms of the new parameters; 

V 2 = [-m\v + X v 3 ] h 



(7.20) 



+ 



3A 



+ [-m 2 1 + \ v 2 ] (G-G + ) + 



X v 2 



(Gl) 



+ \m\\(w 2 w%) + \m 2 3 \(w 3 w£) + \m 2 4 \(w 4 wj 

~\ lV 2 



+ \m 3 \(z 3 z 3 ) + \m A \(zlz 4 ) + 



(44) + 



\m 2 

*„,2 



z 2 z 2 ) 

(z 3 Zi) 



With the value v 2 = one finds the mass spectrum; The SM-like Higgs boson has mass 



2m?, while the masses of the charged scalar bosons are m 2 



\m 



2h 



2 

m ± 



\m 



3b 



2 

m + 

w 4 



\m 



\ I and for the neutral scalar mi 



2 ^3 

4 | aim iui Mic "cuuiai o^aiai — \V(l 2 \. 

Neutrals with equal q can mix, which happens at Ai 7^ 0. Again we describe the 
resulting mass eigenstates by complex fields d and D (m^ < mo), both having q = 3: 



\mH 



|Ai| 2 

± - .„ / (\m 2 3 \ - \m\\) 2 + -^-m^ 



1 



(7.21) 
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where d and D are rotations of z% and z\\ 

d = cos a zs + sin ae~ l13 z%, D = — sin ae l13 z% + cos a z\ 

tan2a = -^- 1 — 9 , mi . — ^ , /3 = argAi (7.22) 
A |mg| — |m|| 

Here we have 

m d < m^i , m^i < m D , m 2 ^ + ± = + (7.23) 

By adjusting free parameters, one can easily make d lighter than w 2 and and 
therefore the lightest among the particles that transform non-trivially under Z 7 . Then, 
the other particles will either eventually decay to d or d* plus SM particles or will be 
stable representing an additional contribution to dark matter. Again, if an asymmetry 
between d or d* exists and if the rate of their annihilation is high, then after the freeze-out 
we are left with the gas predominantly made of d's. 

7.2.2 Higher order terms in the Z-j symmetric potential 

The subsequent microscopic dynamics depend on the interactions between d's and z 2 s 
which follow from the ^-symmetric part of the potential. After eliminating the Gold- 
stone bosons and rewriting the potential in terms of d and D, the interactions between 
d's and z 2 s in the higher order terms in the potential, regardless of the coefficients are 
of the following form 

V 3 = X 1 [hd*d,hd*D,hdD*,hD*D] (7.24) 

+ A 2 [z 2 Z 2 d, Z 2 Z 2 D] + h.C. 

V 4 = X [ddd*d*,ddD*D*,d*d*DD,DDD*D*] 
+Ai [hhd*d, hhd*D, hhdD*, hhD*D] 
+A 2 [hz 2 z 2 d, hz 2 z 2 D] + h.c. 

+A 3 \z 2 w^w^d, z 2 w^w^D, z 2 ddd, z 2 ddD, z^dDD, z 2 DDD] + h.c. 

The relevant terms are hdz 2 z 2 from the X 2 term and dddz 2 from the A3 term. Since 
ma < m Z2 , one- or two-particle processes such as d — > z 2 z 2 , dd — > d*z 2 , and dd — > z 2 z 2 z 2 



are all kinematically forbidden, shown in Figure (7.1): 

Multiple collision kinetics depends on whether m Z2 < 3m^ or not. If z 2 is not too 
heavy, then the "triple semi-annihilation", ddd — > z 2 X$m, is kinematically allowed and 
will create a population of z 2 even if it was absent before. However, z 2 will get de- 
pleted by the semi-annihilation process z 2 z 2 — > d*X$M- So, if one starts with a certain 
concentration of z 2 , d and their antiparticles, then z 2 will die off at a higher rate than 
d's. 

In stationary conditions, the terminal concentrations will be those equilibrating the 
rates of the following 6d tree- level scattering with intermediate z 2 s, 6d — > z 2 z 2 X$m 
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Figure 7.1: Forbidden 1- and 2-particle processes 



d*X' SM and the subsequent annihilation of d*. The net result of this chain will be the 
"7<i-burning process", 7d —> X SM , the bottleneck in this chain being the triple-d process 
ddd — > z 2 X$m, shown in Figure (7.2). 




Figure 7.2: The 7<i-burning processes, with triple-d process bottleneck 



On the other hand, if m Z2 > 3m 



then ddd — > z 2 is kinematically forbidden, while 
the inverse process leads to a quick z 2 decay. In this case, one can still burn d's via the 
tree-level process with intermediate virtual z 2 's, dddd — > d*d*z 2 z 2 — > d*d*d*X$M- The 
net result will be the same 7<i-burning, but the bottleneck process is now the 4d collision, 



whose rate is even more suppressed, shown in Figure (7.3). 
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Chapter 8 
Summary 



In this thesis we reviewed the Brout-Englert-Higgs mechanism of spontaneously sym- 
metry breaking. This mechanism requires the existence of a scalar boson, the so called 
Higgs boson. However, there is no reason to restrict ourselves to only one scalar dou- 
blet, and there are many motivations for introducing more than one doublet. We briefly 
studied the two-Higgs-doublet model in a particular case. 

One of the interesting topics in multi-Higgs-doublet models, is the symmetries that 
one can impose on the scalar potential. It is very helpful to study these symmetries 
in the "orbit space", which was presented here in the general NHDM. We studied the 
particular case of three-Higgs-doublet model orbit space in detail. 

From the analysis of the orbit space we found certain symmetries that are always 
broken after EWSB in multi-Higgs-doublet models with more than two doublets. We 
named them frustrated symmetries because of their resemblance to the phenomenon of 
frustration in condensed matter physics. These symmetries are not specific to doublets. 
They can arise when the representation of the electroweak group has lower dimensionality 
than the horizontal (Higgs family) space, i.e. more than one singlet, more than two 
doublets, more than three triplets, etc. 

Several examples of frustrated symmetries in the three-Higgs-doublet model were 
given; SU(3), octahedral, tetrahedral, and Z3 x Z% symmetry. The case of octahedral 
symmetry was studied further. This symmetry seems to be the largest realizable discrete 
symmetry that can be imposed on the scalar sector in 3HDM. Another interesting feature 
of this case is that it is 2HDM-like, meaning that the mass spectrum of the scalar 
potential with Octahedral symmetry is similar to the one of 2HDM, with a degeneracy 
in the mass of the charged scalar boson, and three different masses for the neutral scalar 
bosons. It would be interesting to see what experimental observables could distinguish 
this model from the actual 2HDM. 

In this thesis we also made a step towards classification of possible symmetries of the 
scalar sector of the NHDM. Namely, we studied which Abelian groups can be realized as 
symmetry groups of the NHDM potential. We proved that they can be either subgroups 
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of the maximal torus, or certain finite Abelian groups which are not subgroups of maximal 
tori. For the subgroups of the maximal torus, we developed an algorithmic strategy 
that gives the full list of possible realizable Abelian symmetries for any given N. We 
illustrated how the strategy works with two small- N examples. We also proved several 
statements concerning the order and structure of finite groups which apply for arbitrary 
N. Finally, we conjectured that, in the NHDM, any finite Abelian group with order 
< 2 N ~ 1 is realizable. 

We extended the strategy of finding Abelian unitary symmetries, to include Abelian 
antiunitary symmetries (with generalized CP transformations) in NHDM. We introduced 
a strategy that could be used for any NHDM, and applied it to the case of 3HDM, where 
the full list of realizable Abelian antiunitary symmetries for 3HDM was found to be 
Z* 2 , Z 2 x Z* 2 , Z 2 x Z 2 x Z|, Z\. 

An obvious direction of future research is to understand phenomenological conse- 
quences of the symmetries found. This includes, in particular, a study of how the sym- 
metries of the potential are broken, and what is the effect of these symmetries on the 
physical Higgs boson spectrum. 

Having found that Z p symmetries are realizable in multi-Higgs-doublet models, we 
showed that these models can naturally accommodate scalar dark matter candidates 
protected by the group Z p . These models do not require any significant fine-tuning 
and can lead to a variety of forms of microscopic dynamics among the dark matter 
candidates (allowing or forbidding semi-annihilation, offering different routes to multi- 
particle annihilation, etc.). 

In particular, we gave explicit examples of dark sectors where the bottleneck process 
leading to depletion of asymmetric dark matter can be a 2-particle, 3-particle or 4-particle 
semi-annihilation. We stress that these models do not require any serious fine-tuning. 
We only ask for the presence of terms invariant under the chosen symmetry group but 
do not constrain coefficients in front of these terms. 

In certain aspects these models resemble the Inert Doublet Model [55j, but in the 
other they rely on symmetry patterns that arise only with several doublets. In this re- 
spect, such models can be viewed as "multi-inert" doublet models although this name 
of course does not completely specify the microscopic dynamics. Exploring the observa- 
tional consequences of each sort of microscopic dynamics is a separate study. 
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Appendix A 
Definitions 



Basic definitions 

• Group. A group G is a set with an operation that combines any two elements 
a and b to form another element. To qualify as a group, G must satisfy four 
requirements known as the group axioms: 

Closure: For all a and b in G, the result of the operation, a b, is also in G. 
Associativity: For all a, b and c in G, (ab)c = a(bc). 

Identity element: There exists an element e in G, such that for every element a in 
G, ea = ae = a. 

Inverse element: For each a in G, there exists an element b = a" 1 in G such that 
ab = ba = e. 

• Group order. The order of a group G, denoted \G\, is the number of its elements. 

• Subgroup. A subset H of G is called a subgroup of G if H also forms a group 
under the group operation. 

• Proper subgroup. If A is a subgroup of B, but not equal to B, meaning there 
exists at least one element of B not contained in A, then A is also a proper (or 
strict) subgroup of B; AC. B. 

• Maximal subgroup. A maximal subgroup H of a group G is a proper subgroup, 
such that no other proper subgroup K contains H. A subgroup of a group is termed 
maximal Abelian subgroup if it is Abelian and is not properly contained in a bigger 
Abelian subgroup. 

• Abelian group. An Abelian group, also called a commutative group, is a 
group in which the result of applying the group operation on two group ele- 
ments does not depend on their order. Example: the groups of even integers 
E = {0, ±2, ±4, ±6, • • • } is Abelian under addition operation. 
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Cyclic group. A cyclic group is a group that can be generated by a single 
element g (the group generator), by applying the group operation on g as many 
times as needed. A cyclic group of finite order n is usually denoted as Z n , and 
its generator g satisfies g n = I, where I is the identity element. Cyclic groups are 
Abelian. 

Unitary group A unitary group of order N, U (N) , is the group oi NxN unitary 
matrices. 

Special unitary group The special unitary group of order N, SU(N), is the 
group of N x N unitary matrices with determinant 1. 

Unitary transformation A unitary transformation is a transformation that 
respects the inner product of two vectors. 

Action 

Group action. A group action is a way of describing symmetries of objects 
using groups. A group G (with operation •) is said to act on a set X when the 
operator o:GxI->I satisfies the following conditions for all elements x G X: 
Associativity: (gi ■ g 2 ) o x = g x o (g 2 o x),Vg 1 ,g 2 e G 

Identity: e o x = x The action of a group on a group is a straightforward gener- 
alization of the action of a group on a set. A group can also act on itself. One of 
the examples of such action is conjugation. 

Conjugation. In a group G, two elements a and b are called conjugate if there ex- 
ists an element g in G with gag' 1 = b. Conjugation preserves the group operation, 
therefore it's a homomorphism. 

Conjugacy class. Conjugation partitions the group into conjugacy classes with 
similar properties. If a is an element of the group G, the conjugacy class of a has 
the form: Cl(a) = gag 1 \ g G G. Every element of the group G belongs to precisely 
one conjugacy class, and the classes Cl(a) and Cl(b) are equal if and only if a and 
b are conjugate. 

Conjugate subgroup. If H is a subgroup of G, with h G H and g a fixed 
element of G which is not a member of H, then the transformation ghg~ l generates 
the subgroup gHg~ l which is conjugate to the subgroup H. 

Homomorphism. A group homomorphism is a map / | G —¥ H between two 
groups such that the group operation is preserved f(gi ■ g 2 ) = f(gi) o f(g 2 ) for all 
gi,g 2 in G, where the product on the left-hand side is the operation in G and on 
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the right-hand side in H. A group homomorphism maps the identity element in G 
to the identity element in H /(e^) = e# . 

Isomorphism. Isomorphism is a homomorphism between two groups with a one 
to one correspondence between group elements. 

Automorphism. A group automorphism / is an isomorphism from a group G to 
itself / | G — > G. f satisfies the following conditions: f(gh) = f(g)f(h) | gi,g 2 G 
G. 

Commutator. The commutator of elements g\ and g 2 of a group G is [<7i,<?2] = 
g\~ x g 2 ~ x g\g 2 - The commutator [<7i, #2] is equal to the identity element e if and only 
if fi'ifi , 2 = <72<7i) that is, if and only if g 1 and g 2 commute. 

Structure 

Center of the group. The center of a group G, denoted Z(G), is the set of 
elements that commute with every element of G: Z{G) = {z G G \ Wg G G, zg = 
gz}. The center is a subgroup of G, which by definition is Abelian. 

Normal subgroup. N is called a normal subgroup of G if for each n G N and 
g G G, the element gng~ l is still inN:A<G <^> V n G A, V g G G , grig' 1 G N. 

Coset. In a group G, with subgroup H, and g an element of G, one defines 
gH = {gh | h G H}as a left coset of H in G, and H g = {hg | h G H} as a right 
coset of H in G. 

Factor (quotient) group. For a group G and its normal subgroup N, the 
quotient group of N in G, G/N, is the set of cosets of N in G; G/N = aN \ a G G. 
The group operation on G/N is the product of subsets, meaning for each aN and 
bN in G/N, the product of aN and bN is (aN)(bN). This operation is closed, 
because (aN)(bN) is a (left) coset itself: (aA^)(6A^) = a(Nb)N = a(bN)N = 
(ab)NN = (ab)N. 

Generating set. A generating set of a group G is a subset S of G, denoted (S), 
such that every element of the group can be expressed as the combination (under 
the group operation) of elements of the subset and their inverses. 

Construction 

Direct product. The direct product of two groups G and H, denoted G x H , 
is defined as follows: The elements of G x H are ordered pairs (g, h), where g G G 
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and h G H. That is, the set of elements of G x H is the Cartesian product of the 
sets G and H. The operation on G x H is defined: (gi, hi) ■ (g 2 , h 2 ) = (<?i-<?2, h 1 -h 2 ). 
Both groups G and H are normal subgroups of the group G x H. 

• Semidirect product. A semidirect product is a generalization of the notion of 
direct product. If N is a normal subgroup of G and H a subgroup of G, we say that 
G is a semidirect product of N and H, G = N x H when G = NH and NnH = e. 

Let Aut(N) be the group of all automorphisms of N, and / be a group homomor- 
phism where / | H — > Aut(N) defined by f(h n ) = hnh},\/h G if, Vn G N. Then 
the operation in G = N x H is determined by the homomorphism /, and is as 
follows; 

* : (N x H) x (N x H) ->■ N x H 
defined by (ni, hi) * (n 2 , h 2 ) = (nif(hi n2 ), h\h 2 ) for ni, n 2 <E N and hi, h 2 G H. 

• Torus. In geometry, a torus is a surface generated by rotating a circle in three 
dimensional space about an axis coplanar with the circle. In topology, a torus is 
homomorphic to the Cartesian product of two circles: Si x Si. An n-dimensional 
torus is the product of n circles. That is: T n = S 1 x S 1 x ■ ■ ■ x S 1 . 

V v ' 

n 

In group theory, a torus in a Lie group G is an Abelian Lie subgroup of G. 
A maximal torus is one which is maximal among such subgroups. Every torus 
is contained in a maximal torus. For example, the unitary group U(N) con- 
tains a maximal torus, the subgroup of all diagonal matrices, of the form: T = 
{diag(e iei , e i02 , . . . , e ie ") : 6j G R}. 
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Appendix B 

Tetrahedral and octahedral 
symmetric potentials 



The potential with tetrahedral symmetry 



The "tetrahedral" potential used in Section 3J2 as an example of a frustrated symmetry 
is defined by the following scalar potential: 



V 



-M r + A r 2 + Ai(r? + r\ + r\) + k 2 {r\ + r\ + r 2 ) + A 3 (r 2 3 + r 2 8 ) 
+A A (r l r 2 - r A r 5 + r 6 r 7 ) 



In terms of doublets, this potential has the form: 



V 



-^(0 t 1 1 +0 t 2 2 + 0^ 3 ) 

Ao 
3 

Ai 
3 

+A 2 



+A 3 
+A 4 



(0^0 + (44> 2 ) + (<f>y 3 ) 

(^ + i0i) 2 + (4<p2) 2 + (4fo) 2 - (0l0i)(40 2 ) - (4<p2)(4fa) - (4<Ps)(4<Pi) 

(Re4<p 2 ) 2 + (Re4<p 3 ) 2 + (Re^f 
(Im0 t 1 2 ) 2 + (Im0 2 3 ) 2 + (Im^i) 5 

(Re0 t 1 2 )(lm0 t 1 2 ) + (Re0 2 3 )(lm0 2 3 ) + (Re<^i)(Im<^i 



This potential is symmetric under the following transformations: independent sign flips 
of the doublets, the cyclic permutation of the three doublets, and their compositions, 
which we show as; 
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Therefore the symmetry group of the potential in terms of these transformation has 24 
elements: 



S 



{e, A B, C, D, E, F, G, G 2 , DA, EA, FA, GA, G 2 A, 
DB, EB, FB, GB, G 2 B, DC, EC, FC, GC, G 2 C} 



The cycle graph of such group is shown in Figure (B.l ), which is the cycle graph of 
a tetrahedron. 



Figure B.l: Cycle graph of the group A4, the symmetry group of a tetrahedron. 



The potential with octahedral symmetry 



The "octahedral" potential used in Section 3.2 as an example of a frustrated symmetry 
is defined by the following scalar potential: 



V 



-M r + A r 2 + A^r? + r\ + r\) + k 2 {r\ + r\ + r 2 ) + A 3 (r 2 3 + r 2 8 ) 



124 



V 



^101 + 0202 + 0303 



Written in terms of doublets, the potential has the following form: 
+Ai L . 

2 



Ao 
3 



+A 2 
As 



(Re^) 2 + (Re0 2 3 ) 2 + (Re^i) 
(Im0i0 2 ) 2 + (Im40 3 ) 2 + (Wj0i) S 



i^O 3 + (050 2 ) 2 + (0j0 3 ) 2 " (0!0l)(0^02) " (0 2 2 ) (0^03) - (tlfcMfa 



a, 



a, 



This potential is symmetric under the following transformations: independent sign flips 
of the doublets, the cyclic permutation of the three doublets, the CP-transformation 
and their compositions, which we show as; 
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Therefore the symmetry group of the potential in terms of these transformation has 48 
elements: 

S = {e, A, B, C, D, E, F, G, G 2 , H, AH, BH, CH, DH, EH, FH, GH, G 2 H, 

DA, EA, FA, GA, G 2 A, DB, EB, FB, GB, G 2 B, DC, EC, FC, GC, G 2 C, 
DAH, EAH, FAH, GAH, G 2 AH, DBH, EBH, FBH, GBH, G 2 BH, 
DCH, ECH, FCH, GCH, G 2 CH} 



The cycle graph of such group is shown in Figure (B.2), which is the cycle graph of 
an octahedron. 
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Figure B.2: Cycle graph of the group S4, the symmetry group of an octahedron. 
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Appendix C 



Detailed study of the symmetric 
minimum (u, v) 



A remarkable phenomenological feature of a 3HDM with Octahedral symmetry is that 
it is 2HDM-like. After symmetry breaking, due to remaining symmetry, this model 
exhibits certain degeneracy in the mass spectrum of the physical Higgs bosons, which 
precisely mimics the typical Higgs spectrum of 2HDM. It would be interesting to see 
what phenomenological consequences distinguish these two models. Here we pick the 
symmetric the minimum point, (v,v,v), and study the potential in more detail. 



Mass matrix eigenstates 

We study the potential at the vicinity if this minimum: 




with 



M 



\/3(Ao + Ai) 
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the mass eigenstates are: 



nt = -^(wt + wt + wt) 



Tin 



-2wJ + wj + wj) 
Pi = ^( h 2 ~ h) 



P2 = ^(^i + h 2 + h 3 ) 

Ps = -^=(-2^ + h 2 + h 3 ) 
V6 

si = ^=(m -vs) 

92 = ^=(?7i +V2 +Vs) 
9s = -j=(-2rii + T] 2 + r] 3 ) 



m 



3 A 2 



< = 



-A lV 2 



2 - 

_ -A1 + A3. 

pi 



-v 



ml = (A + Ai)t» 2 



m P3 = 



-A1+A3 2 
v 



^! = 4(-Ai + A 2 y 



K 3 = 4(- A i+ A 2y 



The three Goldstone bosons that give masses to the and Z are , n 2 and 
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Decay rates 

From the third order terms in the potential we have the decays of the Higgs bosons: 

• p 2 — > n^n^ and p 2 — > n^n^ with vertex coefficient ^(2A — A^v 

• P2 — > P1P1 with vertex coefficient ^(Ao + A 3 )v 

• P2 — > P3P3 with vertex coefficient g (Ai + 2A 3 )t> 

• P2 — > 9i9i and p 2 — > q 3 q 3 with vertex coefficient ^(Ao — + ^A 2 f 

• ps — > nfn^ with vertex coefficient g (Ai — A 3 )t> 

• p 3 — >■ n 3 ng with vertex coefficient — ^(Ax — A 3 )t> 

• p 3 — > q^qi and p 3 — > q 3 q 3 with vertex coefficient — (— A x + A 3 )v 
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• pi — > nfn 3 and pi — > n x n 3 with vertex coefficient ^ (A 3 — A^v 

• pi — > qiq 3 with vertex coefficient ^ (A 3 — A^v 

One could easily calculate that the particles with the same mass p±,p3 and qi,q3, 
have the same decay rate. Following the straightforward method of writing the kinetic 
term, to calculate the couplings to the gauge bosons results in: 

• p 2 -> W + W~ with vertex coefficient ^-g\v 

• p 2 — > ZZ with vertex coefficient ^-(gf + g^) v 

Let's assume A >> A 2 > A 3 > |Ai| > 0. With this choice the mass spectrum 
will have the form m P2 » m qi ^ > m Plt3 > m„ lj3 > m n2 ,m q2 = 0. Therefore the 
decays p 3 — > qi,3qi,3 from the list of interactions above will be forbidden. It appears that 
we would have stable particles in this case, gi j3 and n 13 , which could be dark matter 
candidate. 
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